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PREFACE TO THIRD EDITION 

In preparing this third edition it has not been considered 
necessary to make any change in the lay-out of the text, 
the general form of which seems to have served the needs 
of the physicists and engineers for whom the book was 
written. 

Opportunity has been taken, however,^ to add some 
further notes on vector notation (p. 117), on the more 
formal aspect of the del operator (p. 113) and on Maxwell’s 
equations (p. 118). At the same time a number of minor 
alterations and corrections have been made. I am indebted 
to Professor S. Chapman, F.R.S., and to Professor G. 
Cook, F.R.S., for kindly criticisms and helpful suggestions 
which have enabled me to make the text clearer and more r. 
accurate. For certain historical facts I am indebted to 
the late Dr. John McWhan. 

Glasgow 
Jan, 1945. 


PREFACE 

Vector Analysis is the natural means of expression for the 
three-dimensional problems of physics and engineering, 
because its conciseness and freedom from mathematical 
detail enable the relationships between the various physical 
quantities to be kept clearly in view. Since the pioneer 
work of Gibbs and of Heaviside an increasing number of 
text-books and scientific papers on physical and technical 
subjects have made use of vector methods, until it has now 
become almost essential for any advanced worker in these 
sciences to have some knowledge of vector analysis. Much 
good work can be done with the aid of a very few elementary 
principles. It is the object of this monograph to give an 
introduction to these principles and to explain them from a 
physical standpoint, so that they may be easily available to 
the busy physicist or engineer approaching the subject for 

V 
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the first time. Such workers are usually so much occupied 
by their major task as to lack the time necessary to enable 
them to seek out such principles as they need to use from 
the more comprehensive treatises which aim at mathe¬ 
matical completeness. The monograph is not, therefore, 
intended for the reader with purely mathematical interests, 
whose more rigorous and systematic requirements are fully 
satisfied elsewhere. 

For these reasons the outlook adopted is almost entirely 
physical; geometrical matters and questions of an exclu¬ 
sively mathematical interest are limited to essentials. 
Formal proofs of invariance and conditions of continuity 
in vector processes are replaced by an appeal to physical 
intuition. Purely analytical topics of an advanced kind, 
such as Green’s theorem, are justifiably omitted from such 
an introductory treatment as this monograph aims to pre- 
* sent. The reader who later wishes to amend his knowledge 
in these and other respects will find ample material in the 
standard books listed in the Bibliography. 

The reader’s attention is drawn in particular to two 
features. First, the use of trimetric projection by Gough’s 
method {Engineerings vol. 143, p. 458, 1937) for certain of 
the three-dimensional diagrams, a purpose for which this 
method is admirably suited. Second, the inclusion of a 
chapter giving a brief sketch of the elementary properties 
of tensors and dyadics in their relation to vectors. This is 
a subject which usually puzzles and often repels physics 
and engineering students because of the abstract mathe¬ 
matical way in which it is generally brought to their notice. 

The monograph is based on a course of lectures given 
a few years ago to post-graduate electrical engineering 
students in the Polytechnic Institute of Brooklyn, New 
York. I am particularly grateful to my colleague Dr. A. J, 
Small for his valued assistance in reading the manuscript 
and proofs, and for numerous suggestions. 

Glasgow 
June 1938 
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CHAPTER I 

DEFINITIONS. ELEMENTS OF VECTOR 
ALGEBRA 

1 . Scalar and Vector Quantities. Physical quanti¬ 
ties are divided into two main classes, each with charac¬ 
teristic properties and an appropriate algebra. 

Scalar quantities have magnitude only and do not involve 
direction. Typical scalar quantities are mass, volume, 
density, temperature, electric potential, charge, &c. The 
complete specification of a scalar quantity requires (i) a 
unit of the same kind and (ii) a number stating how many 
times the unit is contained in the quantity. For example, 
to express the mass of a given body we require to know 
whether the unit is the pound, the ton, the gramme, &c., 
and also how many of the chosen units represent the given 
mass. Scalar quantities are manipulated by applying the 
rules of ordinary algebra to their numerical magnitudes ; 
for this reason the algebra of ordinary positive and negative 
numbers is often called scalar algebra and the numbers 
themselves scalars. 

Vector quantities have magnitude and direction. Familiar 
examples are displacement, force, velocity, acceleration, 
stress, electric force, magnetic induction, &c. A vector 
quantity requires for its specification (i) a unit of the same 
kind, disregarding direction, (ii) a number giving the mag¬ 
nitude of the quantity in terms of this unit and (iii) a 
statement of direction. For example, the velocity of a 
moving body is stated by saying (i) that the unit is miles 
per hour, kilometres per second, &c. ; (ii) how many of 
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the chosen units express the magnitude of the velocity 
and (iii) the sense in which the velocity is directed, e.g. 
due north. The combination of conditions (ii) and (iii) 
constitutes the geometric conception of a directed magni¬ 
tude or vector^ quite independent of the kind of vector 
quantity specified by the unit. The directional element 
will prevent the manipulation of vectors by the simple 
numerical algebra applying to scalars ; it is essential, there¬ 
fore, to devise a vector algebra by means of which vectors 
may be handled in a way consistent with the physical 
problems in which vector quantities occur. The laws of 
vector algebra differ in several important respects from 
those of scalar algebra, as will be seen later. 

Vector quantities, though the simplest, are not the only non¬ 
scalar quantities found in physics and engineering. In the theory 
of inhomogeneous strain, in the transformations of space involved 
in the theory of relativity and in the general theory of electro¬ 
magnetic machinery, to mention only a few examples, more com¬ 
plex non-scalars appear. These are the linear vector function, 
the dyadic and the tensor. Each has its own algebra, differing 
radically from the scalar laws and more general in form than 
vector algebra ; indeed, the particular quantity is characterized 
by its algebra. To manipulate these quantities with conciseness 
and ease, higher algebraical methods, multiple algebras and the 
theory of matrices are widely used. In this elementary mono- 
gr./ph it is not possible to do more than sketch the simplest 
properties of these higher non-scalars and their relation to ordinary 
vectors ; the reader will find further developments in the special 
treatises referred to in the Bibliography. 

Scalars are manipulated by six elementary operations, 
namely, addition and subtraction, multiplication and 
division, involution and evolution. In each pair the 
second operation is the inverse of the first. Also, multi¬ 
plication is extended addition, and involution (i.e. raising 
a number to a power) is extended multiplication ; like¬ 
wise, subtraction, division and evolution are successively 
extended inverse operations. All these operations are 
performed by means of certain laws, which will now be 
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stated for addition and multiplication ; they are readily 
extended to the other operations and the reader will find 
formal proofs in any good text-book of algebra: 

(i) The law of association can be written as 

a {h ^ c) {a + i) 

and a {h X c) = {a h) y. c \ 

i.e. the result of adding a to the sum of h and c is the same 
as adding the sum of a and h io c \ with a similar state¬ 
ment for products. In other words, brackets are unneces¬ 
sary in continuous sums and products. 

(ii) The law of commutation which states that additions 
and multiplications can be made in any order ; i.e. 

a + b =- b + a, 
and a y b = b X a, 

(iii) The law of distribution^ 'stating that in compound 
expressions involving multiplication and addition the result 
is equivalent to the sum of individual products taken term 
by term, as in the rules 

ax{b-\-c)^ayh^-axCy 
and (a + b)yc^ayc^hxc. 

The three laws are extended to the inverse operations 
by the rule of signs, and additional rules for the manipula¬ 
tion of indices in involution are also provided. Finally, 
a fourth law, fundamental in the logical discussion of inverse 
operations, is now usually stated, namely, if a x h = 
then cither a = o or b — o. 

While these scalar operations and their laws are familiar 
to all as the working processes of arithmetic and algebra, 
so familiar indeed as to be regarded as self-evident, it is 
quite essential to state them in this formal way since they 
characterize scalar algebra. The essential feature of the 
algebras of vectors, tensors, matrices, &c., is that they all 
violate in some’ way the commutative law for multiplica¬ 
tion and do not necessarily satisfy the condition that 
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either factor in a product is zero if the result is zero. The 
physical meaning of these differences between scalar and 
vector algebra will be explained later. 

2 , Graphical Representation of Vectors. Since a 
vector is the result.of abstracting magnitude and direction 
from a vector quantity, independently of the nature of the 
quantity concerned, we may represent the vector graphi¬ 
cally by a line OA pointing in the direction from O to Aj 
as in Fig. i. The magnitude of the vector is given to a 



Fig. I. —Definition of Equal Vectors 


convenient scale by the length of the line ; the direction 
in space is indicated by an arrow-head marked on the line. 
In the algebra of such vectors it is the custom to distinguish 
them by the use of a distinctive symbol. Thus, the vector 
OA is denoted by the Clarendon letter V ; its scalar magni¬ 
tude is stated by the corresponding italic letter V. A 
vector is unchanged by pure displacement ; it follows, 
therefore, that all parallel vectors of the same length and 
direction are equal; a vector is not localized. Thus in 
Fig. I the vectors V, Vj, Vg, V3 are all equal, a fact stated 
by the relation 

^V^ Vi=.V,^ V3, 

in which the meaning of the sign = is extended to include 
equality of size and similarity of direction. Reversing the 
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arrow-head changes the sign of a vector, such as V4 in 
Fig, I ; then 

indicates equality of size and opposition of direction. 

This definition restricts the meaning of a vector to the repre¬ 
sentation of size and direction only. Consequently, a vector 
cannot represent completely any vector quantity the effect of which 
is changed by simple displacement; many quantities are of this 
kind, most notably force. It is clearly insufficient to specify a 
force by its magnitude and direction alone, i.e. by its vector. 
If the force acts on a rigid body the line of action must also be 
stated, since a change in the line of action alters the torque acting 
on the body. In the case of a deformable body the point of 
action of a force must be specified. These examples are given 
as a warning to the reader not to be tempted to attribute to a 
vector any properties, other than those of size and direction, which 
properly belong to the quantity that the vector is representing to 
scale merely in these particulars. In other words, vector analysis 
must always be used with an eye upon the physical conditions of 
the problem to which it is applied. 

It is geometrically obvious that multiplication of a 
vector by a scalar factor S results in a vector S times the 
size of the first and in the same direction. In particular, 
if V be a vector of unit size in the direction of V, 

V« Fv.(i.i) 

Distinction is often made, particularly in more advanced 
work, between two essentially different classes of quantities 
represented by vectors. The first class includes quantities 
such as force, displacement, velocity, &c., in which the 
vector is drawn in the direction of the quantity concerned ; 
mere linear action in a particular direction is involved and 
the vector is called a polar vector. The second class 
includes quantities such as angular velocity, angular 
acceleration, &c., in which rotary action of some kind takes 
place about an axis. Here the vector is drawn parallel 
to the axis about which the quantity acts ; the length of 
the vector gives the magnitude of the quantity. The 
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direction of the vector is conventionally fixed by the rule 
of the right-handed screw, i.e. when sighting along the 
direction shown by the arrow-head, rotation is taken as 
-clockwise. Such vectors are named axial vectors. The 
distinction is a physical one, there being little differ¬ 
ence in the mathematical treatment of the two classes of 
quantities. Indeed, it is difficult in many cases, particu¬ 
larly in dynamics and electromagnetic theory, to assign 
vector quantities with certainty to classes in this way. 

3 , Addition and Subtraction of Vectors. Consider 
two vectors A and B, shown in Fig, z{a)^ representing for 



Fig. 2.—(a) Vector Addition ; {b) Vector Subtraction ; {c) Vector 
Equation of a Line 


example two successive displacements of a point. Their 
joint effect, sum or resultant is obtained by setting off the 
vector B at the end of A and drawing the vector V joining 
the beginning of A to the end of B. The same result is 
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obtained by starting with B and attaching A to it in a 
similar way, as shown by the dotted lines. Then 

V=^A + B = B+A, . . ^ . (x.z) 

the + sign being understood to mean addition in this 
geometric sense. Thus the sum of two vectors is the 
diagonal of the parallelogram of which the vectors are the 
sides ; such a sum is commutative, i.e. independent of 
which vector is taken first. Subtraction of one vector from 
another immediately follows ; the vector B to be sub¬ 
tracted is reversed and the rule for addition applied, as 
illustrated by Fig, 2(b)y in which V = A-}-(~B) = A—B. 

A useful geometrical application is shown in Fig. 2(c)y 
where A and B are given vectors. If ^ is a variable para¬ 
meter ranging between + co and - oo, the sum of A and 
is the radius vector from O, namely, 

r = A + kB. 

As k varies, the locus of P will be a straight line through A 
parallel to Bj and the expression is the vector form of the 
equation for a straight line. 

The reader must not suppose that all quantities that may be 
represented by a directed line are iiecesmrily vector quantities ; 
the crucial test is whether they follow the parallelogram law of 
addition or not. For example, finite rotation of a rigid body 
about an axis can be represented by a line, of a length propor¬ 
tional to the angle of rotation, drawn in the direction of the axis. 
A second finite rotation about another axis inclined to the first may 
be similarly represented. But their resultant effect cannot be 
found by compounding the two finite rotations by the law for the 
vector sum, as the reader may easily verify in a simple case such 
as the motion of a ball on a plane. Finite rotations are not, 
therefore, capable of representation by vectors ; they i*equire the 
use of tensors and have much more complex laws of manipulation. 
Infinitesimal rotations, angular velocities and accelerations are, 
on the other hand, vector quantities (see p. 46). 

With more than two vectors, such as A, B, C in Fig. 3, 
the sum is obtained by laying off the vectors successively 
2 
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end to end ; then the sum is the closing side V of the 
polygon of which the other sides are A, B and C. But it 
is geometrically obvious that we could start by adding any 



Fig. 3. —Law of Association for Vector Addition 


pair of vectors by the parallelogram rule, and to their 
resultant add the third vector. Hence 
V=A+BtC-(A+B)+C-A+(B+C)=(A+C)+B. . (1.3) 
The sum of any number of vectors is, therefore, associative, 
i.e. the vectors may be added in any desired manner. The 
vectors are not necessarily all in the same plane ; this is 
more clearly shown by Fig. 5. . 

4 . Components of a Vector. By reversing the 
process of addition it is clear that any vector V can be 
decomposed into the sum of n vectors, of which n - 1 are 
arbitrary and the last one closes the polygon. In general 
the vectors are not coplanar and the polygon is not a 
plane figure. 

The most useful instance is the decomposition of a 
vector into component vectors along the three orthogonal 
axes of cartesian co-ordinates. Fig. 4(a) shows the usual 
right-handed system, in which the relative positive direc¬ 
tions along the axes are chosen such that if O'X is turned 
toward O'Y about O'Z through the smaller angle, a right- 
handed screw would advance along the positive direction 
of O’Z ; similar relations hold for the axes O'Y and O'X, 
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always keeping X, Y and Z in cyclic order. Let the origin 
O' be one extremity of a vector V ; draw a rectangular 
parallelepiped with the three edges which meet at O' 


r 




(b) 


Fig, 4. —Cartesian Components of a Vector 

ig along the axes and such that V is the diagonal from 
through the solid figure. If V^c, are the 

tor intercepts along the axes of AT, Y and Z respectively, 
n 


V - V* + + V,. 


. { 1 - 4 ) 
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Now let i, j, k denote vectors of unit magnitude along 
the axes of X, Y and Z respectively; these are the unit 
vectors along the axes and are frequently used in later work. 
Then if Vy^ are the sizes of the vectors Va;, V^, V3, 

from Equation (i.i) ; whence 

V = F,i + .... (i.s) 

In discussing the properties of a vector field we are 
concerned with the magnitude and direction of a vector 
quantity at every point in space. In rectangular co¬ 
ordinates with origin O let a point in space, such as O' in 
Fig. 4(Z>), be specified by the radius vector 

r = + 3;] + ^k. 

Since a vector is not affected by parallel displacement, 
Equation (1.5) is still the relation between the components 
of a vector and the vector itself. Note, however, that as 
we are dealing with a vector field, V is a function of posi¬ 
tion ; that is, V and its components Vy, are func¬ 
tions of the variables x, y, z specifying the position of O' 
relative to O. 

Vector quantities and their vectors are physically inde¬ 
pendent of any system of co-ordinates by means of which 
they may be expressed ; they are said to be invariant. 
Co-ordinate axes are unnecessary, therefore, in vector 
analysis. The resolution of a vector into components, 
given by Equation 1.5, is the connecting link between 
vector notation and that of ordinary co-ordinate geometry 
with which the reader is already familiar. Since a vector has 
three axial components, any vector equation relating vectors 
is the equivalent of three scalar cartesian equations relating 
their components. This conciseness is the great advantage 
of vector methods in three-dimensional problems. The 
fact that vector quantities are expressed directly by a nota¬ 
tion natural to themselves, instead of by a long and artificial 
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f scalar relationships, enables their essential 
Lature to be kept clearly in view. The resolution 
sian components is often useful in the proof of 
used in vector analysis ; in the words of Heavi- 
len in doubt and difficulty, fly to i, j, k advice 
requently follow. 

iple example find the sum of vectors A, B, C. In 
their components, 

= AJ + AJ + AX B - Bji + 

*c = c,i Hh c,j ^ c,k. • 

aiponents along the X axis add directly ; like- 
: along the Y and Z axes. Then, as shown in 


y 



Vyi F^k = A + B + C 

+ (Ay + By + Cy)i 

+ (A^ + B^+C^)k . • (1,6) 

scalar magnitudes in the brackets obey the 
ws of algebra, it follows that vector sums (and 
) obey these laws also. The component of the 
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resultant iri any axial direction is the sum of the individual 
components in that direction ; as is geometrically obvious 
from Fig. 5, 

+ Cz, 

5 . Scalar and Vector Fields. A physical quantity 
can be expressed as a continuous function of the position 
of a point in a region of space ; such a function is called a 
point-function and the region in which it specifies the 
physical quantity is known as a field. Fields are of two 
main kinds, scalar and vector, according to the nature of 
the quantity concerned. 

A typical scalar field, such as the distribution of tem¬ 
perature, density, electric potential or of any other non- 
directed quantity, is represented by a continuous scalar 
function giving the value of the quantity at each point. 
Such a function does not undergo any abrupt change of 
magnitude in passing from any point to another close to 
it, a condition satisfied in all practical cases. The field can 
be mapped graphically by a series of surfaces—such as iso¬ 
thermal, equi-density or equipotential surfaces—^upon each 
of which the scalar has a definite constant value. Such 
surfaces, called equal or level surfaces, are conveniently 
chosen so that in passing from one to the next a constant 
arbitrary difference, is made between the scalars which 
characterize them. It is evident that the level surfaces 
must lie one within the other and cannot cut; for if two 
such surfaces could intersect, the scalar values corre¬ 
sponding to both must hold along their common line, 
which is contrary to our definition. Hence scalar point- 
functions are single valued or uniform at every point. 

A typical vector field, such as the distribution of velocity 
in a fluid or of electric or magnetic field strength, is repre¬ 
sented at every point by a continuous vector function. At 
any given point the function is specified by a vector of 
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definite magnitude and direction, both of which change 
continuously from point to point throughout the field 
region. Starting at any arbitrary place, proceed an in¬ 
finitesimal distance in the direction of the vector at that 
place, arriving at a closely-neighbouring point. Proceed¬ 
ing thence in a similar way, we shall trace out a curved 
line, the tangent to which at any point gives the direction 
of the vector thereat; such a curve is a vector line, line of 
flow ox flux line. To represent the magnitude of the yector, 
at any point on a flux line draw a very small surface per¬ 
pendicular thereto and choose a number of points per unit 
area upon this surface numerically equal to the magnitude 
of the vector. Through each of these points flux lines 
can be drawn. The field is then mapped out by flux lines. 
The direction of the lines is that of the vector function ; 
their density, represented by the number of them crossing 
per unit area perpendicular to their direction, is a measure 
of the magnitude of the vector. It is clear that lines of 
flow cannot intersect, since this would involve indefinite 
direction of the vector at the point where they cut; vector 
point-functions must also be single valued at every point. 

The physical properties of scalar and vector fields will 
be considered in greater detail in Chapters VII and VIII 
respectively. 



CHAPTER II 

THE PRODUCTS OF VECTORS 

1 . General. The ordinary idea of a product in scalar 
algebra, the mere multiplication of a scalar magnitude, 
cannot apply to vectors because of their directional proper¬ 
ties ; nor is it possible to decide by deductive reasoning 
what form the product of two vectors should take. Since 
vectors have their origin in physical problems, definitions 
for the products of vectors must be devised that will be 
consistent with the way in which such products occur in 
applications to physical science. 

As a simple example consider vectors which represent 
a force F and a linear displacement d, their magnitudes 
being F and d and the angle between their directions 0 . 
Products of these two quantities occur in two ways. First, 
the work done by the force is Fd cos 6, a scalar quantity, 
known as the scalar product of F and d. The couple 
exerted by the force has a magnitude Fd sin Q and acts 
about an axis perpendicular to the plane containing F 
and d. It is shown by a vector, called the vector product 
of F and d, drawn along the axis of the couple in a con¬ 
ventional sense to be defined later. We are led, therefore, 
to define two sorts of products, namely, the scalar product 
and the vector product. 

2 . The Scalar Product. The scalar product of two 
vectors A and B is defined as the product of the magnitudes 
of the vectors and the cosine of the angle between their 
directions. As shown by Fig. 6, this is the same as the 
product of the size of one vector with the component of 

H 
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the other in the direction of the first, and this process is 
clearly commutative, since cos d = cos (- 0), i.e. indepen¬ 
dent of the order of the factors/ The scalar product will 



Fig. 6.—Scalar Product of Two Vectors 

he denoted by interposing a dot between the vectors ; 
then 

A*B “ B-A = cos 0 . , . , (^-i) 

Thus, by definition, the scalar product is commutative. 
On p. 3 it has been pointed out that the product of two 
scalars vanishes when either is zero. The scalar product 
of vectors also satisfies this condition, but will vanish in 
addition when the vectors are at right angles. This is the 
first particular in which the laws of vector algebra differ 
from those of scalar algebra. When two vectors are 
perpendicular, therefore, 

AB=o, . (2.2) 

and when they are parallel, 

A-B = ^^.(2.3) 

If B = A, the scalar product of a vector with itself, called 
its self-producty is 

AA = A2 = ^2. ^ ^ ^ ^ (2 4) 

hence to find the size of a vector calculate the square root 
of its self product. 
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Applying these properties to the unit vectors 1, J, k, 
which are mutually perpendicular, 

i-i = j‘k = k i = 0, . . . . (3.5) 

and i2 ^ j2 ^ J5.2 = I. . . . . (2.6) 

Now examine the scalar product of a vector A with the 
sum of two others, B and C, illustrated in Fig. 7(a); the 




(i) 


Fig. 7—(a) Law of Distribution for Scalar Products; (b) Vector 
Equation of a Plane 


three vectors are not necessarily in the same plane. Since 
the resolute of B + C on A is equal to the sum of the 
resolutes of B and C, i.e. OC = OB + BC, 

A-(B + C) = A{OCy= A(OB + BC) = A-B + A-C. ( 2 . 7 ) 
Hence the scalar product is distributive for addition, as 
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is geometrically obvious. By an extension of this process 
it is easy to show that 
(A+B+. . .)-(N+ 0 +. . .) - 

A-N+A-O-i-. . . B-N+B‘ 0 +. . . +.(2.8) 

A particular case is of geometrical interest. In Fig. 7(6) 
let p be the perpendicular from O upon a plane which 
passes through the extremity of a given vector c. If r 
is the radius vector to any other point in the plane, then 
it is obvious that the resolutes of c and r upon p must 
be equal, i.e. 

p-r = p-c or p'(r - c) = o 

is the equation bf a plane perpendicular to p through the 
end of c. From Equation 2.2, (r — c) must lie in the plane 
normally to p, as is obvious. 

Referring again to Fig. 2(a), the scalar product enables 
the usual formula for the resultant of two vectors to be 
readily found. Their sum is 

V-A + B. 

Take the self product on both sides ; then from Equations 
2.4 and, 2.8, 

v-v - V2 = (A + B)-(A + B)= A-A + 2 A-B + B-B, 

or -h zAB cos 6 + 

The rectangular components of a vector are found at 
once by taking its scalar product with the unit vectors, as 
is evident from Fig. 4 ; then 

VI = V,, VJ = F, and V-k = F,. . . (2-9) 

The scalar product assumes an important form in 
rectangular co-ordinates. Writing from Equation 1.5, 
AJ + Ayj + Ajis. and B ^ + BJ + 

Equations 2.5, 2,6 and 2.8 give 

A-B = {A^i + AJ + + Byi + B,k) • 

^ +AyB^ +AJB^. , (2.10) 
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That is, the scalar product of two vectors is the sum of 
the products of their components along each of the co¬ 
ordinate axes. If, for example, A is a force and B is a 
displacement, the total work done is equal to the sum of the 
works done by the components of force and the correspond¬ 
ing components of displacement, as is physically obvious. 

3 . Line and Surface Integrals as Scalar Products. 
In Fig. 8, let OA be any curve drawn in a vector field and 



dl an element of length along it at any point P. Let V 
denote the vector at P in a direction making an angle 0 with 
that of the length element. Then V‘dl is the product of 
the length of the element and the component of V in its 
direction, i.e. tangentially to the curve ; thus 

V-dl=Vcos 6 dI. 

If V varies in magnitude and direction from point to point 
along the curve, the integral 

cos edl , . 

Jo Jo 

is defined as the Ime integral of V along the curve OA, 
Such integrals are of very frequent occurrence. For 
example, if V is a force and di an element of the path of a 
particle along any curve, the line integral denotes the work 
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done in displacing the particle from O to A. Again, if V 
is the electric field strength, i.e. the force on unit charge in 
an electric field, then the line integral expresses the 
potential difference between the points O and A. Further, 
if V is the velocity at any point in a fluid and the integral 
be taken round a closed curve, the integral is called the 
circulation of the fluid. In an electromagnetic field the 
line integral of the electric force round a closed path is the 
electromotive force in the path ; in a purely electrostatic 
field the e.m.f. is zero. 

Consider now any element of area ds upon a surface 
drawn in a vector field, and let V be the value of the vector 
aft the middle of the element. Draw the positive* normal 
of unit length n on the element and let 0 be the angle 
between n and V, as in Fig. 9. Then the component of 



V perpendicular to is V*n F cos 0 and the flux of V 
through the element is V-n ds. The integral of this taken 

* If the surface is closed, n is drawn outward from the enclosed 
volume. If the surface is unclosed, then n is drawn always on the 
same side of the surface, the normal on the opposite side being 
negative. 
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over the surface is called the total flux or surface integral of 
V through the whole surface i.e. 

V cos d ds. , . (2.12) 

The physical meaning of such surface integrals is very 
simple. Suppose V to denote the vector velocity of a 
moving fluid in which a fixed surface 5 is drawn. At any 
point V'li ds denotes the amount of fluid passing normally 
through the surface element ds in unit time, and this com¬ 
ponent alone must be considered since the tangential com¬ 
ponent of V necessarily contributes nothing to the flow 
through the surface ^element. If the flow through an ele¬ 
ment is in the sense of the positive or outward normal 
it is counted as positive ; if in the opposite sense, as on 
the element ds^, it is negative. The integral of these 
normal contributions expresses the aggregate flow of fluid 
through the whole surface in unit time. If the surface is 
closed, positive total flux diverges from the enclosed volume 
while negative flux converges upon it. Should the amount 
of flux entering by some elements be balanced by that 
leaving by other elements in such a way that the total flux 
is zero, then either there are no sources or sinks of fluid 
within the enclosed volume or their sum is zero. Similar 
ideas apply to other fluxes, e.g. of electric or magnetic 
induction, of heat, &c. 

4 . The Vector Product. The vector product of two 
vectors A and B is defined as a vector having a magnitude 
equal to the product of the magnitudes of the factors and 
the sine of the angle between their directions, and a direc¬ 
tion perpendicular to the plane containing A and B. The 
sense of the product along this perpendicular is defined by 
the right-hand screw rule ; if the vector A is turned to¬ 
wards B through the smaller angle, the necessary rotation 
must be clockwise when sighting along the positive normal 
unit vector n in Fig. 10. That is, the rotation needed to 
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move A to the position of B and the positive direction of n 
are related in the same way as the rotation and translation 
of a right-handed screw ; this direction of n is that 
of the vector product. It follows from this definition that 
a change in the order of the factors in a vector product 
reverses the sign of the product, since sin (- 0) = ~ sin 0. 

AxB 



BxA=-AxB 

Fig. io. —Vector Product of Two Vectors 

Hence the vector product of two vectors is not commutative 
and the order of the terms must be strictly maintained, 
a further difference between vector and scalar algebras. 
The vector product will be denoted by a cross between the 
factors ; then 

AxB« — B X A <= AB sin 0 n, . . (2-13) 

where n is a positive normal of unit length drawn from the 
plane containing A and B. The product vanishes not 
only when either factor is 2:ero but also when the com¬ 
ponent vectors are parallel, again differing from scalar 
algebra. When two vectors are parallel, 

A X B = o, * . . . (2.14) 

and when they are perpendicular, 

A X B = n, .... (2.15) 
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in which case the two vectors and their product are mutually 
at right angles. Applying these ideas in particular to the 
cartesian unit vectors i, j, k, ^ 

ixi = ixj=kxk-o, . , . (2.16) 

and 

ixj=-jxi=k, jxk—kxj=i, kxi=-~ixk-j. . (2.17) 

The strict cyclic order in these important results should be 
observed. 

5 , Vector Area. The magnitude of a vector product 
is AB sin 6, i.e. the area of the parallelogram with sides 
^ and B and included angle d. Since the shape of the 
figure is not specified by the definition of the product, any 
plane area of amount AB sin 0 with its positive unit normal 
n can be taken to represent a vector product. By an 
extension of this idea we are led to the notion of vector 
area. A plane area, such as s in Fig. 11(<2), can be regarded 
as possessing both magnitude and direction. Its magnitude 
is the amount of the area and its direction is that of the 
normal to its plane. The sign to be attributed to a vector 
area is defined with reference to the order in which it is 
traced out as viewed from an external point such as O. 
If the direction of tracing is counterclockwise as seen from 
this viewpoint, the positive direction of the vector area is 
along the unit normal n, the directions of tracing and of n 
being related by the right-hand screw rule. Then s = sn 
is the definition of a vector area s. Vector areas can be 
resolved or added just like other vector quantities. 

In Fig. ii(Z)) a tetrahedron is shown with vectors 
S2, S3, S4 drawn to represent the vector areas of its four 
triangular faces, the outward normal being regarded as 
positive. Resolve these areas upon any plane, e.g. one of 
the faces of the solid. Then some of the projections will 
be positive and others negative, the sum of them all being 
zero. Hence the total vector area of. a tetrahedron is 
zero. This result follows at once from physical consider 
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tions, if the tetrahedron is supposed to be drawn within a 
fluid which is in equilibrium under hydrostatic pressure. 
Each face experiences a force normal to its plane and 
r^roportional to its area. Since the fluid within the solid 
figure is in equilibrium with that outside it, the resultant 
of the forces on its faces is zero ; hence also is the sum of 
tifie vector areas, since the pressure is the same on all faces. 
TT'his hydrostatic demonstration also applies to any form 
of solid figure and thus generalizes the theorem. Geo- 
ixietrically, any polyhedral figure may be divided up into 



{a) . {b) 

Fig. II.— (a) Vector Area ; (/j) Vector Area of a Closed Surface 


tetrahedra ; every surface introduced into the interior of 
tHe polyhedron appears twice, once with a positive and 
once with a negative normal. Hence for any polyhedral 
sxarface the total vector area is zero. By making the faces 
vanishingly small and indefinitely increasing their number 
w^e approach a closed curved surfiice over which 

jjis = o. 

6. Application to Vector Products. Comider now 
trlie vector product of a vector A with the sum of two 
ot:Hers, B and C. In Fig. 12 draw a triangular prism with 
its parallel edges in the direction of A and its end faces as 
tr-iangles with sides B, C and B + C. The vector areas 
of the triangular end faces are JB x C and x B, which 
3 
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cancel; the remaining vector areas are B x A, C x A 
and A X (B + C) in the sense of the outward m'rmals, 
their sum being zero. Thus 
Ax(B + C) + BxA + CxA = o 

i.e. Ax(B+C)—BxA-CxA=AxB+AxC . . (a.i8) 



Fig. 12.—Law of Distribution for Vector Product 


Hence the vector product is distributive, but the order of 
the factors must be strictly observed. By repeating the 
process it is easy to show that 

(A+B+. . .)x(N+ 0 +. . .)= 

AxN+AxO+. . . BxN+BxO+. . (2.19) 

for any number of vectors. 

These rules enable the vector product of two vectors to 
be expressed in rectangular components. Using the nota¬ 
tion of Equation (1.5), 
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on applying the rules for vector products of unit vectors 
and the distributive law, Equations 2.16, 2.17 and 2.19. 
This result is more easily remembered by writing it in the 
form of a determinant ; thus, 


A X B - 


i i k 

Ay A, 

By B^ 


(2.20 bis) 


7 . Products of Three Vector^. The vector product 
of two vectors B and C being a vector, can give both a 
scalar and a vector product with a third vector A. There 
are, therefore, two triple products, namely A-(B x C) and 
A X (B X C), which occur frequently in physical applica¬ 
tions. 

The scalar triple product has a simple interpretation 
illustrated by Fig, 13, (B x C) is a vector normal to the 


B><C 



plane of B and C, with a magnitude equal to the area of 
the shaded parallelogram. The scalar product of A with 
(B X C) is the product of this vector area and the projec¬ 
tion of A along (B x C) ; i.e., A-(B x C) is the volume of 
the parallelepiped which has A, B, C for its edges. Any 
face of this solid can be taken as the base ; hence three 
equivalent expressions for the volume are 

A-(B X C) « B-(C X A) « C-(A x B), . (2.21) 
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cyclic order of the factors being maintained to retain the 
volume with positive sign.* Since the order of terms 
in a scalar product is immaterial, Equation 2.1, these rela¬ 
tions are equivalent to 

(A X B)-C = (B X C)-A = (C X A)-B ; . (2.21) 
so that the dot and cross may be interchanged at will. 


y 



Fig. 14. —^The Vector Triple Product 


Three vectors, therefore, have six identical scalar triple 
products which may be written concisely as 

i j k A^AyA, 

[ABC] = {A^i+Ayi^A^yiy By = B^ By B^ . (2.21 b. 

Cy Ci Cj; Cy Cj 

When three vectors lie in a plane, the volume of the 

* Scalar quantities of this kind, in which the sign depends on 
the proper cyclic arrangement of the component vector factors, 
are called by mathematicians * pseudoscalars to distinguish them 
from true scalars, which do not change sign when the reference 
axes are changed from a right-handed to a left-handed system. 
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parallelepiped is zero ; hence the condition for vectors to 
be coplanar is that their scalar triple product vanishes. 

The vector triple product A x (B x C) can also be illus¬ 
trated geometrically, as in Fig. 14. The vector (B x C) 
is normal to the plane containing B and C. Likewise, the 
vector A X (B X C) is normal to the plane containing 
A and (B x C), i.e. in the same plane as B and C. Take 
the X axis along B, the Y axis at right angles thereto in the 
plane of B and C, and the Z axis along (B x C). Then 
in terms of their axial components the three vectors are, 
by Equation 1,5, 

A-= A^i + AJ + B = BJl and C = + Cyi, 

since this choice of axes makes = Cg = o. Ex¬ 

panding the vector product (B x C) by Equation 2,20, 

(B X C) = ; 

and also by the same rule, 

A X (B X C) = AyB^Gyi - A^B^Cyl 
But from Equation (a. 10) 

A*C = Aj^C^ + AyCy and A-B = AJB^y 

so that 

A X (B X C) - AyCyB - AACyi 

= {AyGy + A^Cg;-)^ — Ag-Gg^BJl — A^Bg^Gy} 

= (Aa + - AAiC^^ + Cyi). 

Substituting from above, 

A X (B X C) = (A-C)B - (A-B)C. . . (2.23) 
Each term of tl\e product involves the external factor A 
in a scalar product, first with the extreme and then with 
the middle factor. Hence the value of a vector triple 
product is entirely different if the order of the factors be 
interchanged. 

Products of more than three vectors do not often occur 
in physical applications and when encountered they are 
easily reduced by use of the preceding theorems. The 
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results serve to confirm what is already shown by Equations 
(2.21), that vectors do not follow the law of association 
for multiplication, p. 3. 

8. Summary. Chapters I and II give all the essential 
rules for the manipulation of vector algebra. By extending 
the meaning of the terms addition and subtraction to the 
geometrical processes explained on p. 6, sums and differ¬ 
ences of vectors follow the ordinary laws of commutation 
and association, as shown by Equations 1.2 and 1.3. It is, 
therefore, immaterial what order or manner of grouping 
of the vectors is followed when adding or subtracting them. 
Vectors and scalars are so far in agreement, but important 
differences appear when products are considered. 

The scalar product of two vectors is commutative and 
distributive, Equations 2.1 and 2.7, and vanishes when 
either factor is zero, exactly as in ordinary algebra. It 
differs, however, in vanishing when the vectors are per¬ 
pendicular. The vector product of two vectors differs 
considerably from the ordinary rules. It is not commuta¬ 
tive (Equation 2.13), it vanishes when the vectors are 
parallel, and it is distributive only if the order of the 
factors is maintained unchanged (Equation 2.18). For 
more than two vectors the non-commutative nature of the 
vector product causes the ordinary law of association to be 
not obeyed, as has been demonstrated in Section 7 of this 
chapter. 

Mathematically, vectors can be defined as quantities 
that require for their manipulation a non-commutative 
algebra which has the agreements with, and differences 
from, ordinary scalar algebra that have been summarized 
above. Vector algebra is only one of many non-commu¬ 
tative algebras known to mathematicians, but it is among 
those that have the greatest practical utility. 
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THE DIFFERENTIATION OF VECTORS 


1 . Scalar Differentiation. Let V be a vector which 
js undergoing a continuous change of magnitude anci 
^Jirection. In Fig, 15 ( 5 V is a small increment which gives 


d^r 6v 
w w 



Frc. 15.—Scaliir Differentiation of a Vector 


tzlie new value V (!)V for the vector. Suppose tliat V 
is a vector function of a scalar variable t ; then wlien t 
eli,?inges from t to t h- rV, V becomes V r^V. The ratio 
3 ^/^t is the average rate of change of V with /, and as 
3 t becomes vanishingly small the ratio attains a liinuing 
value which is the rate of increase of V, i.e, 


(ly 

dt 



us 6 t becomes zero. 
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This is the derivative of V with respect to the scalar 
variable t. As V varies, its extremity moves over a curve 
and the derivative is a vector in the direction of the tangent 
to this curve at each point. Second, third and higher 
derivatives are obtained by similar arguments, by analogy 
with the ordinary ideas of scalar calculus. 

When V is expressed in rectangular co-ordinates, 

V=L,i+F,j + F,k; 

are now functions of t. Since i, j, k are 
constant vectors, it follows that 


£V 

dt 


dV^ 

dt 


I + 


dV^,. 

dt 


■J + 



• • ( 3 -i) 


and similarly for higher derivatives. 

A most important physical instance is when V is a 
vector of displacement and t is time ; the extremity of V 
is then the path of a moving body. The derivative 
dW/dt is the velocity along the path at any instant, and 
d^'^/dt^ is the corresponding acceleration. 

2 . Differentiation of Sum6 and Products. If V 
is the sum of two vectors A and B, both of which are 
functions of then a change of ^ to t dt makes 

' V + < 5 V * (A + 5 A) + (B + < 5 B), 
so that 5 V = ( 5 A + dB. 

Dividing by 6 t and proceeding to the limit gives 




dt dt' 


• . • ( 3 - 2 ) 


so that the operation of differentiation is distributive, as 
in ordinary calculus. 

For the scalar product V = A-B the increment in t gives 

V+(5V=(A+6A)-(B+aB)=A-B+6A-B+A-5B+(5A-dB, 

expanding the scalar product by Equation 3.8. Subtract* 
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ing A'B, dividing by dt and proceeding to the limit, 
neglecting second-order terms, 

^ / 4 -TIN -r. 4 , X 

^'li' ■ ' ■ 

The order of the factors can be changed at will, i.e. the 
operation is commutative. If B = A we have the useful 
form for the self product, 

= = = 3 ^. ( 3 . 3 «) 

For the particular case of a constant vector =: == constant, 

defining the position of a point on the surface of a sphere of radius 
A, dA/dt is zero and Equation 3.3^2 leads to A‘(dA/dt) ^ o. 
Hence dA/dt is perpendicular to A. If a point moves on the 
surface of a sphere its velocity is always normal to the radius 
vector, as is physically obvious. 

For the vector product V = A x B expansion by Equa¬ 
tion 2.19 gives 

V + ( 5 V = (A + ( 5 A) X (B + ( 5 B) - A X B + dA X B 

+ A X ^B + ( 5 A X 5 B, 

leading to 

Z(A’<B)-f xB+Axf, . . ( 34 ) 

•in which the order of the factors must be strictly maintained. 

Putting p ^d/dt the reader can easily verify that for 
triple products 

^[A-(B X C)] - ^ A-(B X C) + A-(>B x C) + A-(B x pQ), 
and 

^[A X (B X C)] X (B X C) + A X (j 5 )B X C) 

+ A X (B X ^C). 

In general, differentiation in vector analysis is seen to 
follow the .same rules as in ordinary differential calculus, 
except so far as vector algebra differs from scalar algebra 
in the non-commutative property of a vector product. 
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3 . Partial Differentiation. These simple properties 
of differentiation as applied to vectors can be extended to 
partial derivatives when a vector is a function of more 
th^ one scalar independent variable. The most useful 
case is that of a vector V which is a function of the car¬ 
tesian co-ordinates Xy z of s. point in space. If 3^ and z 
remain constant while x increases, the partial derivative 
dV /dx denotes the rate of increase of V with respect to x. 
Likewise changing y and z alone gives the partial deriva¬ 
tives BY/dy and BY/Bz, denoting the rates of increase with 
respect to y and z respectively. If now Xy y and z change 
simultaneously by differential increments dxy dyy dz, the 
total change or total differential of V will be 


dV = ^dx + ^dy + —ds 
Bx By Bz 


( 3 - 5 ) 


which is of frequent .occurrence in physical applications of 
vector analysis. If r = jtfi + + ;sk is the radius vector 

from the origin, then its differential increment is 


dr = dxi + dyj + dzk. 


Equation 3.5 may be symbolically written as 

If we now define the operator V by 


^ ^dx'^hy 



it is easy to verify that the scalar product of V > regarded 
as a kind of vector, with dr gives the operator in square 
brackets. Thus, 

dY = (^/•dr)Y .(3.5 bis) 

The operator V is of immense importance in physical 
applications of vector analysis, wherein it -appears in 
association with both scalar and vector operands. These 
uses will be fully explained in Chapters IV and V. 
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THE OPERATOR V AND ITS USES 

1. The Operator V. The differential operator V was 
introduced by Sir William Rowan Hamilton and developed 
by P. G. Tait; it is of central importance in all three- 
dimensional physical problems. The symbol V was origin¬ 
ally named ‘ nabla ’ after a harp-like ancient Assyrian 
musical instrument of similar shape ; other writers have 
called it ‘ atled \ i.e. ‘ delta * reversed. It is now usual 
to adopt the term ‘ del ’ introduced by J. Willard Gibbs. 
In cartesian notation 

del = V=i^ + j^ + k^, . . . ( 4 - 1 ) 

which may be applied as a directive differentiator either 
to a scalar or to a vector function of space. Again, treat¬ 
ing'the differentiators in V as scalars, we may formally 
regard V as a vector which can have either a scalar or a 
vector product with other vectors. In vector analysis 
there are three fundamental operations with V which are 
of physical interest. If *5 is a scalar function and V a 
vector function of space, these operations are (i) V*S', where 

V acts as a differentiator ; (ii) V'V, and (iii) V x V, where 

V is treated as a formal vector. (See also p. 113). 

2 . The Gradient of a Scalar Field. On p. 12 it 
has been shown that certain physical quantities, such as 
temperature or electric potential or any such non-directed 
quantity, can be represented from point to point in space 
by a scalar point-function, iS, of the co-ordinates. The 
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entire scalar field can be mapped out by level surfaces, 
upon each of which the scalar function S has a constant 
value. Consider two such surfaces very close together 
and examine a small portion of them in the neighbourhood 
of a given point A on the surface characterized by the 
constant value S of the scalar function ; the second surface 
is specified by a constant value S + dS, This is shown, 

un/'t 

norma! 

n 



Fig. i6. —Gradient of a Scalar Point Function 

much magnified, by Fig. 16. If r is the radius vector 
from the origin to the point A, any point such as B in the 
second surface is given by r + dv. The least distance 
between the surfaces will be AC^ in the direction of the 
unit normal vector n 2Lt A and of length dn. 

If dr be the length of AB, the magnitude of the rate of 
increase at A oi S in the direction of AB will be dS/Br 
when the two surfaces are vanishingly close together. 
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This rate of increase becomes‘greatest in the direction of 
the unit normal n, i.e. along AC, when it has the value 
dS/dn, Note that 


ds 

dr 


ds . 

= 7;^ cos 0, 
on 


Hence, if n is the unit vector normal to the level surface 
at any point in a scalar field, the vector ndS/dn gives the 
greatest rate of increase of S at the point in magnitude 
and direction." This vector is called the gradient of >5 at 
the point and is written 

gradS = ^n.(4.2) 


Thus, the gradient of a scalar field is a vector field, the 
vector at any point having a magnitude equal to the most 
rapid rate of increase of S at the point and in the direction 
of this fastest rate of increase, i.e. perpendicular to the 
level surface at the point. 

A simple physical example will fix the reader’s ideas. 
Suppose S is the potential in an electric field due to static 
charges. Then the electric force at any point is in the 
direction of the greatest rate of decrease of potential, i.e. 
normal to the equipotentiahsurfaces, and has a magnitude 
equal to that rate of decrease. That is, the electric force 
is — grad aS* and is obviously a vector field. 

It is clear from the way that the idea of gradient has 
been introduced that it is an intrinsic property of a scalar 
field and is, therefore, a physical notion entirely inde¬ 
pendent of any particular system of co-ordinate axes. 
In other words, the operation denoted by grad is invariant, 

2 a. The Operation VS, Consider now the vector 
represented by 


_ ^ dS. dS. dS, 
“ dx^ dy^ dz^ 


in rectangular co-ordinates. The vector rates of increase 
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of S in. the directions of the axes of x, y and z are idS/Bxy 
idS/dy and VidS/dz, Their sum will be a vector with 
the magnitude and direction of the total or most rapid 
rate of increase of *S. It remains to show that this expres¬ 
sion is equivalent to the gradient as just defined. To do 
this, take the scalar product on both sides of the gradient 
equation with an element of radius vector dr^ i.e. resolve 
dr in the direction of the normal at A in Fig. 16. Then 

(grad Sydr = ^^n-dr = ^dr cos d = ^dn = dS, 


since dS/dn is the total normal rate of change of S, Now 
in rectangular co-ordinates, 


05 , 05 . 05 , 

dS ^ ^dx + ; 

ox ay dz 


hence, 

(grad sydr = + ^^dz 

r JO r-7 O d5. 05. 05., . v 

so that grad 5 = V s = + _j + -k, . . {4.3) 

showing that the operations grad and V applied to a 
scalar point-function are identical. Also since the gradient 
in a scalar field is independent of the nature of co-ordinates, 
so also is the operator V ; it is a mere convenience, there¬ 
fore, to express it in terms of x, y and z. We may regard 
V as a directional differentiator which, when applied to a 
scalar function of space, gives the greatest rate of change 
of the function in magnitude and direction at every point, 
i.e. derives from a scalar field its vector field of gradient. 

We have here an important case in which a vector field 
is derived from a scalar field by the process of finding 
the gradient of the latter. It does not necessarily follow 
conversely that all vector fields can be expressed as the 
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gradient of a scalar function. Let be a vector which 
is derived from a scalar S in the form 

- grad S - 

In Fig. 17 draw any path, such as that marked Path r, 


V^=ffradS 



Fig. 17. —Line Integrals in a Lamellar Field 

between two points A and B in the vector field and let 
make an angle 0 with the element dl of the path ; 
then the product of the length of the element and the 
component of in its direction is, from p. 18, 

cos Odl= V,,v/1. 

The path is traced out by the extremity of a radius vector 
r from the origin ; then dl is equivalent to dr and we may 
write, from p. 36, 

Y^'dl « - (grad Sydr = dS. 

Hence the line integral of V^i from ^ to 5 is 



VECTOR ANALYSIS 


3S 

where Sj^ and Sjj are the values of the scalar at the ex¬ 
tremities of the path. Since only the end values are 
concerned, the line integral has the same value inde¬ 
pendently of the path followed. 

Now take a closed path, consisting of Path i from 
A to B followed by Path 2 from B to A. From Equation 
(44), since the limits of integration are reversed for the 
second path, the value of the line integral along it is 
> hence round a closed path 

^(grad = 0 .... (4.5) 

Summarizing, when a vector field can he expressed as 
* the gradient of a scalar fields the line integral of the vector 
taken between two points is independent of the path 
followed and is equal to the difference between the values 
of the scalar at its ends ; further, the line integral round 
any closed path in -such a vector field is zero. 

A vector field derived from a scalar by the relation 
== grad S is sometimes called a scalar potential fields 
S being the potential of V^. Since space is divided up 
into layers or laminae by the level or equipotential surfaces 
of the function is also called a lamellar vector. 

Again, since the essential property of such a vector is 
that its line integral round any closed path is zero, the 
vector field is called a non-curl fields a term which wil 
be explained in Section 4. 

A simple example is the electrostatic field of chargee 
conductors, which can be specified by a system of equi 
potential surfaces. The vector field of electric force i 
derived from the scalar potential as its negative gradient 
i.e. the electric force is in the direction of the greatest rat 
of fall of potential and has a magnitude equal to that rate 
The line integral of the electric force between two point 
is the potential difference between them and is independer 
of the path taken ; it represents the work done in moviu 
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a unit charge of positive electricity from one point to the 
other. Round a closed path the work done is zero, ex¬ 
pressed by the statement that the electromotive force 
vanishes round any closed path in a static electric field. 
The reader will be able to construct other physical examples 
of lamellar vector fields, e.g. in the flow of heat with tem¬ 
perature isothermals or in gravitational attraction with 
level surfaces. 

3 . The Divergence of a Vector Field. In Fig. i8 



Fig. i8. —Divergence of a Vector Point Function 


let V be the value of a vector function at the middle of 
an infinitesimal element of volume with sides dXy dy and 
dz parallel to the axes of y and z. The vector V has 
axial components of magnitude Vy and V^, To fix 
ideas the reader may think of V as a vector giving the 
velocity of a moving fluid in magnitude and direction. 

Consider the two fiices of the volume element, each 
with area dy dz perpendicular to the axis of X, On the 
left-hand face the value of the x component of the vector 
at the middle of the face becomes 

dV, 
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and may be taken as the value all over the face when this 
becomes vanishingly small. Similarly, on the right-hand 
face the x component is 






+ 


We now define the flux through any face as the scalar 
product of the vector area of the face and the vector V ; 
i.e. as the product of the area of the face and the normal 
component of the vector upon it, see p. 19. The flux 
is positive when the component of V and the outward- 
drawn normal on the face of the element are in the same 
sense. Then the excess of flux leaving the element over 
that entering it in the X direction is 


dz=^^dx dy dz. 

In the hydromechanical case this represents the net volume 
of fluid passing per second in the X direction. By similar 
reasoning the contributions parallel to Y and Z are 

dy dz and dy dz. 

The total net flux diverging from or leaving the element 
is, by Equation 2.12, 


dx dy 



The amount of this flux per unit volume is defined as the 
dwergmce of the vector V and is written 


dVx dV, 
div V ^ + 


dV, 


Since the divergence is the atnount of flux, it is essentially 
scalar. 

If the divergence exists at a point in a fluid, whether 
liquid or gas, and is positive, it expresses the rate at which 
fluid is flowing away from the point per unit volume thereat. 
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Hence, either the fluid is expanding and its density at the 
point is falling with time, or the point is a source at which 
fluid is entering the field. When the divergence is nega¬ 
tive it gives the rate at which fluid is flowing toward the 
point per unit volume. In this case either the fluid is 
contracting and its density rising at the point, or the 
point is a negative source, a sink^ at which fluid is leaving 
the field. Since most practical liquids are almost incom¬ 
pressible, the existence of divergence in them means the 
presence of a volume distribution of sources or sinks rather 
than changes of density. In the case of non-material 
fluxes, such as those of the thermal, electric or magnetic 
fields, the existence of divergence means the presence of 
a source or sink of flux at the point. For example, in the 
electric field, positive, divergence means that there is 
positive electricity at the point; in the thermal field the 
point is either a source of heat or a place where the tem¬ 
perature is falling. 

When the divergence is everywhere zero, the flux enter¬ 
ing any element of space is exactly balanced by that leaving 
it and we may write, 

div V = o,.(4.7) 

which is true in many practical problems. In a fluid 
this means that there can be no sources or sinks in the 
field, nor can its density be changing ; i.e. the fluid is 
incompressible. If the fluxes entering and leaving an 
element are equal, none can have been generated within 
it ; the lines of flow of the vector V must either form 
closed curves (cf. the magnetic field of a current), or 
terminate upon bounding surfaces (cf. the electric field 
in a condenser), or extend to infinity. A vector which 
satisfies this condition is called solenoidal (from a Greek 
word meaning a tube). 

3 fl. The Operation V‘V. Consider now the scalar 
product of the operator V and the vector V. Expressed 



42 


VECTOR ANALYSIS 


in rectangular co-ordinates the rules for scalar products, 
Equations 2.5, 2.6 and 2,8 give 




dx By dz 


that is, 


VV-div Y 


(4.8) 


Now the operation V been shown on physical grounds 
to be invariant, i.e. independent of any system .of axes. 
Hence the idea of divergence is also invariant. This is, 
indeed, physically obvious, since the amount of flux enter¬ 
ing or leaving per unit volume at a point is clearly a con¬ 
ception quite independent of any system of co-ordinates 
that may be used to express the shape and position of 
the volume element. 

4 . The Curl of a Vector Field. It has been shown 
on p. 38 that when a vector field can be derived as the 
gradient of a scalar field, the line integral of the vector 
taken round any closed path is zero ; this result is true 
no matter what size or shape the path may have. A 
vector field satisfying this condition is known as a lamellar 
field, and is of a special, though very important, class. 
Many vector fields occur in physical problems, however, 
in which the closed-path line integral is not zero and 
which cannot, therefore, be expressed as the gradient of 
a.scalar point function; it is to an important property 
of these more general fields that we now give attention. 

Consider a very small region of such a vector field, 
several lines of flow in which are shown by Fig. ; 

the portion is chosen small enough for the lines to be 
regarded as nearly straight and parallel. Into this field 
put a small plane area, shown for convenience as a rect¬ 
angle. 'When the area is perpendicular to the field, as 
in position i, none of the field is directed along any part 
of the bounding edge of the area ; the line integral round 
it is zero. In position 2, with the area parallel to the 
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field, since the value of the vector along the upper edge 
•is assumed to be different from that along the lower edge, 
the line integral round the boundary has a finite value. 
Similar argunients apply to intermediate positions ; the 
value of the line integral depends, therefore, upon the 
direction of the normal to the area relative to the field. 



Fig. 19.—Cud of a Vector Point Function 


i.e. upon the orientation of the given small vector area at 
the region considered. 

In general, if we put a small vector area of any shape 
at any point in a vector field and compute the line integral 
of the vector V around its bounding edge there will be 
an orientation of the area for which the line integral is 
greatest. The amount of this maximum line integral ex¬ 
pressed per unit area is called the curl of the vector field 
at the point, and is given the vectorial sense of the positive 
normal drawn on the small exploring area when in the 
position giving this greatest integral. Some writers use 
the term rotation (briefly rot,) to denote the same con¬ 
ception, since the curl is associated in hydromechanics 
with the rotation or spiii of a fluid in which the particles 
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have angular velocity. Physical examples illustrating the 
meaning of the curl of a vector field will be given later ; 
the meaning depends very much upon the nature of the 
problem. 

To calculate the curl in terms of its cartesian com¬ 
ponents, take three infinitesimal rectangular areas inter¬ 
secting mutually at right angles at a point where the 
vector field V has components of magnitude Vy, V^, 
as in Fig. 19(6). Taking the positive normals to the areas 
along the positive directions of the X, Y and Z axes 
respectively, the circular arrows indicate the positive senses 
in which their boundaries must be traversed to accord 
with the right-hand screw rule for vector areas. 

As an example, take one of these areas such as abed in 
Fig. I9(^) with sides dx, dy, its normal being along the 
axis of Z. Since the rectangle is very small, the numerical 
value of the component of V at the middle of any side may 
reasonably be taken as the average value along that side ; 
the arrows show the directions in which the components 
act. Since Vy, are functions of the co-ordinates 
(^j y, z) of the middle of the rectangle, the average values 
along the four sides ab, be, dc, ad will be respectively - 


-1 

^ X 


dv. 


.dV, 


dV^ 


,dV„ 


and Vy -h^dx. 
.Around the contour abed the line integral is, therefore, 


+ [^^ 2 , + - (Vy - \-^d:^dy, 

that is, 



Since the area of the element is dx dy, the bracketed term 
is the magnitude of the component curl of the vector field 
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taken about the Z axis. Giving it the sense of unit vector 
k, therefore, we may write 


curl-j V = 




By following a precisely similar method with the two 
remaining rectangles we find the component curls about 
the y and X axes to be 


Wi ^ 

dy dzj 


curly and curly V 

Adding the three components gives 


dz 


dx y dx dy 


( 4 - 9 ) 


which can be conveniently written in the form of a deter¬ 
minant as 

i j k 

curl V = d/dx d/dy djdz . . (4.9 hii) 

y. Vy Vz 

4 a. The Operation V x V. Take the vector product 

V X V of the operator V and the vector V, expressing 
the result in rectangular co-ordinates. Using Equation 
Z,20 for expanding a vector product, 

V ><v= 




dy 


i|; + k£) x(F.i + F,j + F,k) 
dz J y dz dx j ) ' 


that is, 

V X V = curl V ... (4.10J 

The operation is, of course, independent of the system 
of axes. 

5 . Simple Examples of Curl/ The meaning of the 
operation curl may now be illustrated by a few simple 
examples. First consider a rigid body which is rotating 
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with an angular velocity co about an axis OA^ 0 being a 
fixed point in the body. Then any point P, Fig. 20, 



Fig. 20. —^Curl and Angular Velocity 
moves in a circular path about OA with a tangential linear 
velocity cor sin 6, where r is the distance from O to P, 
in vector notation the angular velocity is a vector co 
drawn along the axis in a sense related to the rotation by 
the right-hand screw rule ; the co-ordinate of P is the 
radius vector r. Then the tangential velocity at P is 
CO X r, which has magnitude cor sin 6 and is perpendicular 
to the plane containing co and r.* If in addition the 

* To verify that it is physically correct to treat an angular velocity 
as a vector it is only needful to show that two such velocities 
compound as vectors, see p. 7. Let cOi be the angular velocity 
about an axis through 0 ; then the tangential velocity at P is 
X r. Now let to 2 be the angular velocity about a second 
axis through O inclined to the first ; then P has a tangential velocity 
due to this rotation of 0)2 x r. But linear velocities compound 
vectorially, hence the total linear velocity at P is 

toi X r + toa X r = (cOi + 0)2) x r = to x r, 
by the distributive law, Equation 2.18. Hence the motion is the 
same as that due to an angular velocity equal to the vector sum of 
Ae components. A similar theorem is obviously true of 
infimtesimal angular rotations and of angular accelerations ; it 
iS not true of finite angular rotations-, ^ 
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whole body is moving with a linear velocity Vq in any 
direction, the total velocity at P is 

V == Vo + to X r, 

and curl V = V xV = V xVq + V xtoxr. 

Now since Vq is constant for all points in the body and is 
thus independent of x, Zy V x Vo is clearly zero. 

The angular velocity to is also a constant vector for all 
points and can be written as 

to = co^i + co^j + co^k, 

where its components co^ are independent, there¬ 

fore, of the coordinate r of the point P. Writing 

r xl + y] + ark, 

and using Equation 3.20 makes 

to X r = {o)yZ - co^y)i + (co^x - C0j^z)i + (co^y - a)yX)k, 

Again, using Equation 2.20 with Equation 4.1 and remem¬ 
bering that since co^, are not functions Xy 3/, z their 

derivatives are zero, it is easy to verify that 

V X (to X r) = 2C0^i + 2COyj + 2co^k « 2t0. 

Finally, therefore, 

curl V = 2to 

and to = J curl V.( 4 -ij) 

Hence when a rigid body is in motion the curl of its 
linear velocity at any point gives twice its angular velocity 
in magnitude and direction. (See also p. 113.) 

Now examine the motion of a fluid and consider what 
may happen to an infinitesimal volume within it. The 
volume can have three kinds of motion simultaneously : 
(i) It may be moving with a linear velocity of translation 
as a whole, (ii) If the fluid is deformable the element 
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may charge its shape, (iii) It may be in rotation. At 
any instant the little volume may be regarded as a rigid 
body ; the curl of the velocity of the fluid at the point 
where the element is situated gives twice its angular 
velocity, the instantaneous axis of rotation being that of 
the curl. The nature of the rotational motion will be 
made clearer by Fig. 21 (<2), which shows two positions of 




Fig. 21. —Rotational {a) and Irrotational {h) Motion 

a small portion of the fluid in movement about an axis 
at O. It is clear that the portion has rotated ; and if 
%yery elementary portion of the volume immediately round 
1 ^ has rotated by the same amount, curl V would give 
twice the angular velocity of rotation about O, By con¬ 
trast, Fig. 2i{h) shows a smafl portion of fluid which, in 
its motion about O, does not rotate ; hence there is no 
curl V and its angular velocity is zero. 

If a motion is such that the velocity has a curl, the 
ultimate particles of the body are in rotation with an 
instantaneous angular velocity. The motion is described 
as rotational or vortical. If, on the other hand, the par¬ 
ticles do not rotate, there is no curl and the motion is 
irrotational or non-vortical. 

Another example of an entirely different kind occurs in 
the magnetic field of a conductor carrying a steady current. 
At any point in the field put a very small plane area and 
turn it into such a position that the line integral of the| 
magnetic force taken round its boundary is the greatest 
possible ; this value expressed per unit area is the vector 
curl H, i.e. the magneto-motive force per unit area at 
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the point. If the point is within the material of the 
conductor at a place where the vector of current density 
is c, then this will be the total current passing normally 
per unit area through the elementary path when the line 
integral round it is greatest. In electromagnetic units, 
curl H = 4?rc, so that the curl of the magnetic force at 
any point is proportional to the current density thereat 
and has the same direction. For a point in the field 
external to the conductor there is no current density and 
curl H = 0. 


6. Divergence of a Vector Product. In certain physical 
problems, e.g. in calculating the flux of energy in an electromag¬ 
netic field, it is required to find the vector product of two vectors 
and then to work out the divergence of the result. Treating this 
at present as an exercise in manipulation of vector operations, we 
are to find div (A x B). Using Equation 2.20, write 

V «A +{AxBy—AyBx)^ 
^Vxi+Vy} + Vzk. 


Now apply Equation 4.6 and 



Using Equations 2.10 and 4.9 gives 


div (A X B) = B- curl A — A* curl B . . (4.12) 

7 . Divergence and Curl of ^S^A. In later work we require 
the divergence and curl of a vector field which is itself the product 
of a scalar field S and a vector field A at every point. The com¬ 
ponents of the product are SAx^ SAy and SA^,, 

Using Equation 4.6 to find the divergence, 

div ■— + ^{SAy) + -^{SA^ 




: 5 div A + {AJ. + Ay] + Aiiky(^i + 
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which from Equation 4.3 is 

div SA = S div A + A* grad aS , , . (4*13) 

Using Equation 4.9 bis the curl can be written 

1 j k 

curl SA = d/dx d/dy djdz , 

SA^ SAy SAz ' 

the X component of which has the magnitude 



Using Equation 4.9, the first bracket is the x component of curl A ; 
the second bracket is the x component of (grad S) x A, as may be 
verified from Equations 4.3 and 2.20. Working out the other 
components gives finally 


curl SA t= ^ curl A + (grad S) x A 

s S curl A — A X grad S . , (4.14) 



CHAPTER V 


FURTHER APPLICATIONS OF THE OPERATOR V 


1. The Operator*-div grad. In Equations 4.2 and 
4.3 it has been shown that if < 5 * is a scalar function of 
position in space, then grad xS* == V * 5 ' is a vector function 
having the magnitude and direction of the greatest rate of 
increase of S ; this vector is the gradient of S and the 
flux lines of grad S cut normally through the level surfaces 
of S. Since grad 5 is a vector it can have a divergence ; 
using Equations 4.8 and 2.4 this scalar quantity is 

div grad xS - V •(V*S') - V‘VC^S') - \ 7 ^S. 

From Equations 4.1 and 4.3, 


div grad S 



.3 ,d\/dS. dS, 



Expanding by the usual rules for scalar products, Equations 
2.5, 2.6 and 2.8, 


div grad S => 


^ ^ ^ 
dx^ dy^ 


• • ( 5 -i) 


and the operation 

02 02 02 

div grad = + = . . (5-2) 

which is known as Laplace's operator. Examples of its 
use will be found in later chapters. 

2 . The Operator curl grad. Since grad xS is a vector, 
it is possible to calculate its curl and thus interpret the 
operation curl grad applied to a scalar point-function. 
Now from Equation 4.9 his and above, since the corn¬ 
el 
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ponents of the vector for which the curl is to be found 
are dS/dx, dS/dy and dS/d%, 

i j k 

curl grad S d/dx d/dy d/dz = zero . (5.3) 

'=Vx(V^)“ dS/dx dS/dy dS/dz 

This result follows at once from physical considerations. 
On p. 38 it has been shown that = grad iS is a lamellar 
vector, of which S is the scalar potential. The charac¬ 
teristic feature of a lamellar field is that the line integral 
of round any closed path is zero. Since the curl of 
a vector field is a particular kind of closed-path line 
integral, curl = curl grad S is necessarily zero as 
is lamellar. A typical non-curl field of this class is the 
electric force due to static electric charges. 

3 . The Operator with Vector Operand. A 
vector point-function V may be expressed in terms of its 
rectangular components V^^}, V^k in the usual way. 
Since Kc, Vy and are each scalar functions of position 

yy ^)y Equation 5.2 applies to them all. Then if 

V = + v,i + V,K 

V^v = -f S/Wyi + . . (54) 

which is a result of great importance in electromagnetic 
theory and hydromechanics. 

4 . The Operator grad div. If V is a vector field, 
div V is a scalar field, which, therefore, has a gradient. 
This new vector, of which div V is the potential, is neces¬ 
sarily lamellar because curl grad on scalar operand is 
zero ; it may be written as 


grad div V=V (V •V)= 


.a .a , a 


dV.^. dVy 
dx dy dz 


y a^^ dxdy dxdzj 





dy^ dydz 


\dxdz dydz dz^ j 


(5-5) 
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in cartesians. For example, if V is a fluid velocity, div V 
gives the scalar rate at which the density at a point is 
changing per second ; grad div V gives the magnitude 
and direction in space of the greatest rate of space increase 
of the density. In an electric field V is the electric force, 
div V the density of space charge and grad div V the 
greatest rate of increase of charge at any point. 

5 . The Operator div curl. If V is any vector field, 
curl V is also a vector field, for which we may calculate 
the divergence and the curl. Dealing with the first of 
these, using Equations 4.6 and 4.9, 

dy 


div curl V =« V '(V V) 


dz ) 


dyy dz dx 




= zero 


(5.6) 


dzy dx dy 

Hence the curl of any vector field is itself a solenoidal 
field, i.e. one in which the lines of flow form closed curves, 
since it has no divergence (see p. 41). 

6. The Operator curl curl. The y and z components 
of curl V are 


[ dz 


dx / 


and 


/dVy_ 

\dx 


dy 


k; 


and hence by Equation 4.9 the x component of curl curl V 
will be 

curL curl V 

Expanding the factors, adding and subtracting dW^/dx^ 
gives 


"1® 

dVA 

_ l(dV^ 


dy\ dx 

dy) 

dz 

dx ) 


curL curl V 


( 


d^ ^ d^ 

dx'"^ dxdy'^dxds. 


; 


dx^'^dy^^dz-^^ 


F.i, 


which from Equations 5.5 and 5.4 may be identified as 
the X components of grad div V and V^V respectively. 


Permuting i, j, k ; x, y, 


and Vy, in cyclic 
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order gives curl^ curl V and curl^ curl V respectively ; 
adding the three components 

curl curl Y = V x (V x Y) == grad div V - V^V, (5.7) 
so that the operation 

curl curl - grad div - • • . (5.8) 

a theorem which is much used in various problems. 

7 - The Classification of Vector Fields. We have 
distinguished so far between only two kinds of vector 
fields. The first is characterized by zero curl and can 
be expressed as the gradient of a scalar potential, p. 38 ; 
it is, therefore, known as a lamellar or scalar potential 
field. The second has zero divergence ; there are no 
sources of vector flux in the field, the flux lines are closed 
curves and the field is termed solenoidaly p. 41. Since 
the operation div curl on a vector operand yields zero, 
Equation 5.6, a solenoidal vector can always be expressed 
as the curl of another vector ; see below. 

In practice there are four kinds of vector fields, having 
different associations of curl and divergence ; they may 
be classified as follows : 

(i) curl Y =« o and div V = o ; Since the field is irrota- 

tional or lamellar the first condition makes V = grad S, 
The second condition makes the field also solenoidal; 
from Equation 5.1, div Y = div grad o, which 

is Laplace’s equation. Typical fields of this sort are the 
electric force in free space due to static charges on boundary 
conductors, or the irrotational motion of an incompressible 
fluid. See Chapter VII. 

(ii) curl V ~ o but div V o : Here the field is still 
irrotational but is not solenoidal. Hence V o, which 
is Poisson’s equation. The field must contain a spatial 
distribution of sources of flux. Examples are the electric 
field of a volume distribution of charges, e.g. electrons in 
a thermionic tube, or the gravitational force inside a mass. 
See Chapter VII. 
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(iii) curl V =*= o but div V = o. The field is rotational 
and cannot have a scalar potential, but it is solenoidal. 
Let A be a vector such that V == curl A; then from 
Equation 5.6, div V = div curl A — o as required. Equation 
5.7 makes curl V = curl curl A = grad div A — V^A. The 
nature of the vector A is a matter of choice ; suppose, 
therefore, that we desire it to be solenoidal, then div A = o 
and 

curl V = - V^A ^ o. 

This is similar in form to Poisson’s equation but contains 
a vector instead of a scalar operand. For this reason A 
is termed the vector potential of V. Important typical 
cases are the magnetic field within a conductor carrying 
a steady current and the rotational motion of an incom¬ 
pressible fluid. See Chapter VIII. 

(iv) curl V 4= o and div V 4 = o. This is the most 
general type of vector field possible, found in the rotational 
motion of compressible fluid and in some branches of 
electromagnetism. The field can be expressed as the sum 
of two fields ; the first is a lamellar vector, which has no 
curl but can have divergence, and the second is a solenoidal 
vector, which has no divergence but can have curl. Writing 

V = grad S + curl A, 
div V = div grad 5 == V ^* 5 ^ 4 = o, 
from Equations 5.1. and 5.6. This is Poisson’s equation 
and determines S. Using Equations 5.3 and 5.8 and 
choosing A to be solenoidal, as before, 

curl V = curl curl A ^ - V 4 = o, 

which determines A. Thus a general vector field may be 
decomposed into two fields, one being lamellar, with 
scalar potential and the other solenoidal, with vector 
potential A. This is sometimes known as Helmholtz's 
theorem. 

Fig. 22 illustrates diagrammatically the characteristics of vector 
fields near a point. Supposing the vectors to illustrate the 

5 
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motion of fluid particles, Fig. 22(a) shows neither divergence nor 
curl. Fig. 22(h) shows divergence alone ; fluid is streaming from 
the point, at which the density is falling ; alternatively the point 
is a source. Fig, 22(c) shows fluid in rotation about the point; 
there is curl but no divergence, and the fluid is incompressible. 
Finally, Fig. 22(d) illustrates particles in a compressible fluid in 
rotation ; there is both curl and divergence. 


t 

- 4 

Q 

divergence 

no divergence divergence 

cur! 

no cur! 

only 

only 

and curl 

(«) 

(h) 

w 
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Fig. 22 ,—^Typical Characteristics of Vector Fields 


8^ Two Useful Divergence Formulae. In the applications 
of vector analysis two special forms of the divergence formula, 
Equation 4.13, are often used. 

Let 'U and W be two scalar point-functions, grad W being a 
lamellar vector field with W as its potential. Then in Equation 
4.13 put U for and grad W for A, giving 

div (U grad IF) = UV W + V W 

= UV (grad t 7 )-(grad W) . (5.9) 

Again, let S he a scalar field and A = curl B a vector field, 
B being the vector potential of A ; then Equation 4.13 becomes 

div (S curl B) « ^ div curl B + (curl B)*(grad S). 

But Equation 5.6 makes the first term zero ; hence 

div (S curl B) = (curl B)-(grad S), . . (5.10) 

9 . The Vector Field grad (k/r). One of the most 
important lamellar fields is that which has a scalar potential 
S ^ k/ty k being a numeric and r the length of the radius 
vector T = oci + y] + from the origin to any point P 
in the field. The level surfaces are concentric spheres 
with a common centre at the origin ; the potential increases 
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toward infinity as the origin is approached. From Equa¬ 
tion 4.3, 

grad S “ grad (k/r) = V(V'') 

*= — k(xi + jyj + <3'k)/(A:^ + 

Simplifying, 

grad (k/r) = - kr/r^ = - kr^/r^, . . (5.11) 

where ri - r/r is unit vector in the direction of r. As 


spherical 



shown by the minus sign the gradient is directed radially 
inwards, as drawn in Fig. 23 ; the lines of flow of the vector 
field grad {kjr) are thus radial lines cutting perpendicularly 
through the spherical level surfaces and converging upon 
the origin. Excluding for the moment the origin itself, 
at which h/r and grad (A;/r) both become infinite, calculate 
the divergence at any other point by means of Equation 5.1. 
Then the reader can easily verify by direct partial differen¬ 
tiation of k/r that 

SJ^k/r) = div grad {k/r) *=0, . . (5.12) 

so that the potential satisfies Laplace’s equation, p. 54, 
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everywhere except the origin, wlicre the divergence is 
infinitely great. 

Hence when the scalar potential is i/r, its gradient is 
--ri/r®. This is the familiar case of the unit point source 
or sink common in all branches of tlicoretical physics. A 
potential of this kind is known ns a Netotmiati potential^ 
since it is that of a point unit gravitating mass at the origin ; 
the gradient gives the attractive force per unit moss placed 
at P, this force varying with the inverse square of the 
distance. Many other phenomena also follow a similar 
law. A negative electric charge at the origin attracts n 
positive cliarge at P; if both are of unit magnitude, 
Coulomb’s law shows Equation 5.11 to be true. Again, 
in hydromechanics the velocity of the fluid toward a point 
sink—where fluid is abstracted from the field—varies in¬ 
versely as the square of the distance and satisfies the same 
equation. 

On p. 56, if unit mass attracts a similar mass at unit 
distance with unit force, ^ » i. Following the convention 
of p. 13, if a single flux line per unit areij of tlie level 
surface is drawn perpendicular thereto to represent unit 
value of the gradient at unit distance, it follows that the 
total flux of grad (I/r) is —471, since unit sphere has a 
surface of 471 and the flux is inward. If the quantity of 
attracting matter at O is ft, the flux is 47z:ft. I’he subject 
will be further discussed in Chapter VII. 



CHAPTER VI 

THE THEOREMS OF GAUSS AND STOKES 

1 . The Divergence Theorem of Gauss, Consider 
a <ilo3ed surface s drawn in a vector field V ; the normal 
flux tlirougli an element of vector area =« nds^ ii being 
the outward drawn normal, is V-rfe ■= V*n^f^, as explained 
on p. 19. The surface integral of such elements will 
give the total norrrial flux through the surface in Hie form 
given by Equation 2.12, 

V= II 

Referring to Fig. 24, an elementary volume dv within 


V 


Fio. 24.—Gauss*8 Theorem 

the surface s is shown ; a small cube has been taken for 
conycniencc, but any shape will suflicc, 'Phe size of the 
element is much exaggerated, merely for ease in lira wing 
the diagram. The total flux diverging from this volume 
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18 div Vdv, V being the vector at the mitltllc of die element, 
and this Hux Jiaa been coiiiputcd on ji. 40 by calculating 
the aOrface integral of tlio normal componcjus of V 
through the Biirfaccs bouiuling the volume, bor the face 
abed the positive direction of the A? component of V and 
the outward-drawn normal arc in the am no sense and the 
flux is positive. A contiguous cube, ahowu dotted, Iipa 
the same component of V acting through the coininon 
face—since V and its componoiils arc iisBUincd to bo 
continuous, as also are the derivntives of these qiiaiililics 
but its outward normal, shown dotted, is in the opposite 
sense aJid the flux ie negative* The siirfiicc fluxes tluougli 
the common face cancel. Tiiis argument can be apjflictl 
by adding further elements along die X axis until wo reach 
an clenicnt which has one faco in die aurfneo shown 
shaded; this alone makes n contribution to the aorjiial 
flux through j. Applying the snmo trcatmciu to the 
surfaces of all elements of volume throughout the space 

enclosed by j> wo arrive at a total flux V*rJs, but 

flt the same time wo have integrated div Vdv throughout 
the enclosed volume aiul this also incnaiircs the total flux. 
Equating them, 

j'|v-is JJv-Wi- |j^|div JJ 


which is Gam*s theorem of divergence. 

Although the deduction of this theorem m not n rigid innlhc- 
imtical proof—many such can ho found in the larger icxt-hookK— 
It ia based on eelf-ovidcnt phyaical facta. CojiRldcr, for exam pie, 
0 closed aurfdco dra^vn within n fluid wlilcii ia niovinu at n given 
point with vector velocity V. Tho total nmoiiiit of fluid jxinaing 
per second through the surface can be found in two c(|,iiividont 
ways. First, by calculating V-i/a, l.c. tho product of an eJejuent 
of Borfacc nnd the component of velocity pcr|ictu[jcular to it, for 
every element of tho surface and ncldiriR all ihu coEitrihijtionH. 
Second, by investigating the divergence of a volniiio element, i.o* 
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tho exccyjB of eoiircea of fluid over the sinks contained in it per 
unit volume, and integrating cliv Vdv througli the enclosed 
volume. The two results nro phyaicnlly equivalent, einco tho 
excess of fluid leaving the surface over that entering it must bo duo 
to the net amount of fluid injcctctl into tho portion of tho field 
within I by the aggregoto of Boiircca and sinks. 

2 . Gauss*8 Thoorom and tho lavorso Square Law. 
In Fig. 2^{a) let (fe be an elementary portion of a aurfiice, 
n being unit positive normal upon it. From a point O 
draw a conical pencil touching the boundary of r/s anti 
let Ti be unit vector in the direction of the radius vector 
r « rFj from O to the element. Draw splieres of radii r 
and unity respectively; If do) is the area cut by tho 
cone from the surface of tho sphere of unit radius, 
cos O/r^ ; doo is called the solid angle sub¬ 
tended by ds2it O. The solid angle ie regarded as positive 
when the angle 0 is acute and negative when obtuse. 
The total solid angle subtended by a spherical siirfaco 

IS 

Fig. 25(^f) shows any closed non-splicrical surface. 
When O is inside tho eiirface an element at i, 2, or 3 
contributes a positive solid angle. A cone joining O to 
elements 4, 5, and 6 shows’tlint these subtend equal Bolid 
angles at O, but while the contributions of 4 and 6 arc 
positive that of 5 is negative. In goncrnl, therefore, any 
small cone from O cuts the surface an odd number of 
times, equivalent to a single contribution of doi, IIciicc 
the entire surface s subtends at O a total solid angle of 
47r, When tho vertex of the cone ia at ‘O' outside the 
surface, the surface is cut an even number of timcH ^viL]^ 
alternately negative and positive equal cojitributicms, ns 
shown at 7, 8, 9, 10. Honce the anrface s aubteiulH a 
solid angle of zero at an extcrjuil point. 'These two facls 
arc a purely geometric consequence of tlie definition of 
solid angle. 

Consider now the Newtonian potential field of a p(}int 
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unit mflss or other similar physical entity placed either 
at O or at O'. Taking its potential aa i/r, Equation 5.11 
shows that we can put in Equation 6.1 

V - grad {i/r) - - r^/r^ 

and V*^s « - {y^/r^yi\ds - - C03 Qds/r'^ - - da)y 

noting that Vi and n are unit vectors with an angle 6 
between them. Then 



Hence the total closed siirfaco integral of the vector 
grad (i/r)) e.g. of the total flux of force from a point unit 
mass, is “ 4JC wheji the mass is inside tlic surface and 
zero when outside it. The mmus sign means that the 
flux is inward, see also Fig. 23. In both cases the flux 
is independent of the position of the mnsa, which may bo 
anywhere inside or outside the surface. This form of 
Gauss's theorem 19 cxlrcincly useful aJtd will be used in 
Chapter VII. 

3 . Stokes’s Theorem. In n vector field V draw any 
unclosed surface or cap having for its hounding edge n 
given closed curvo> and calculate the lino integral of V 
round the closed curve when it is traced out counter¬ 
clockwise as shown, by tlio boiuKlnry nrrows in Fig. 26. 
Then from Equation 2.11 the value of this integral will bo 

the circle indicating that a closed path hns been travoracd. 

By means of a system of intersecting linos drawn in 
the unclosed surface let it be divided up into infinitesimal 
surface elements. Consider the shaded clement of vector 
arca nr/^, n being the positive unit nonual upon it. 'The 
boundary of the element is-traced out comitci-clockwise, 
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this sense and the direction of n following the ususl right- 
hand screw rule, see p. 22. At the middle of ds the 
vector field will have ft curl and the scalar product of this 
with the vector area will give the component of the line 
integral of V round the boundary of ds, l.e, n*curl Vds, 
as explained on p* 43. A similar process can be applied 
to all surface elements, tracing them all out in the same 
way. Since V and its derivatives, and hence curl V, are 



Fio. 26,—Stokes’s Theorem 


all supposed to be finite and continuous functions, the 
line integral along the common sides of contiguous elements 
will cancel, leaving only a gontribution from sides that lie 
in the given closed curve. Adding the contributions of nil 
elements yields the closed line integral, but in obtaining 
it a surface integral of n-curl Vds has been simultaneously 
made ; hence, 

"" I Jn-curl yds, . . • (6-3) 

This is Stokes's theorem and states that the tangential lino 
integral of a vector function round any closed curve is 
equal to the normal surface integral of the curl of that 
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function over any unclosed surface which has die curve 
for its bounding edge ; rigid mathematical proofs taking 
Account of all the necessary continuity conditions would 
he out of place here and will be found in the usual text¬ 
books. The converse theorem is necessarily true, i.e. if 
the line integral of V round a closed curve is equal to the 
surface integral of A over any open surface bounded by 
the curve, then A curl V, 

"When curl V is zero, Equation 6.3 shows that the closed- 
path line integral of V vanishes ; this has been shown on 
p. 38 to be a characteristic feature of a lamellar field. 

4 , Invariance of Divergence mid Curl. Gnusa^a injicl 
Stokcfl’a tlieoreina may be regarded from eevcral points of view. 
Their physical meaning haa already been strcaacd, the first theorem 
providing alternative ways of expressing the flux of a vector through 
n closed surface, the second giving equivalent Ideas for the line 
integral of a vector round a closed curve, Both results are 
necessarily scalnr quantities. From another standpoint wo may 
look upon them as useful onalyticnl transformations which enable 
certain integrals to he more simply expressed. Equation 6.1 
shows that the volume integral of a vector function can he ex- 
pressed as a surface integral over n closed surface confining the 
volume, Equation 6.3 reduces the integral of a vector over an 
open surface to one taken round its boundary edge. 

Returning to Equation 6.1, confine the volume integral to the 
space within an infinitesimal clement rfu and the surface integral 
to its enclosing surface, then 

The divergence at a point is the value of tide surface integral per 
unit volume as the clement is made vanishingly small, i.e. 

VrfsJ as dv becomes zero. 

Similarly, apply Stokes's theorem to an infinitesimal piano area 
fit a given point, turning the area until the line integral of V round 
its boundary has the largest value ; then n and curl V are in the 


div V ^ Lim, 
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Brtine direction and their sctdai product is the inagnitudo of ih^ 
curl, say curl V, giving 

(b Vv/I - curl V ds, 

J ma£. 

Hence, as ds is inade vanishingly small, 
w/V-Lim. 

. and has the direction of n in the position of mnximiim line Integral. 
The ideas of divergence and curl have been shown on physical 
grounds to be conceptionR independent of any syatem of co¬ 
ordinates in which they may be expressed ; they arc iitvariatd. 
The formulae derived above can be regarded os a purely mathe¬ 
matical juBtifiention of this property, since no co-ordinates ore 
involved in their deduction. Some mathematical treatises, indeed, 
treat these formulae as,definitions of div and curl, building up the 
theory of vector analysis thereon. 



CHAPTER VII 

THE SCALAR POTENTIAL FIELD 

1, General Properties. Let iS be a scalar point-func¬ 
tion which may be mapped out in space by n scries of 
level surfacesy upon each of which the scalar has a definite 
but different constant value. These surfaces divide up 
the region of space into a series of layers or laminae. 
Associated therewith is a vector field directed every¬ 
where normal to the level surfaces, i.e. in the direction 
of the greatest rate of increase of S at any point and having 
a magnitude equal to that rate of increase. This is ex¬ 
pressed by Equation 4.3, 

V ,9 - grnd 5 « V*?. 

Fig, 27 illustrates tlicsc conditions for small portions of 


grad S 



Fia. 27.—The Scalar Potential Field 


two infinitcsiinaily-close level surfaces, characterized by 
values of S and S -v dS respectively. 
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Referring again to Fig. 17 and Equation 4.4, it has been 
shown that such a vector field hna the essential property 
that the value of the tangential line integral of Y.s 
any path joining two points A and B is the same and equal 
to the difference between the values o( S nt B and A, i*e, 

- s^. 

Hence round any closed path the line integral of is 
zero ; in particular round a path of infinitesimal size the 
integral also vanishes, i.e, curl is zero. 

Because of the relation between and V5 is called 
a lamellar or scalar potential vector fields S being potent iaL 
Since the curl is zero, the field is also called a non-curl 
or irroiational field, 

2 . The Inverse Square Law. Point Sources. In 
practice the greatest interest is found in lamellar fields 
based upon the law of inverse squares. Consider a point 
source of vector flux, around which the level surfaces are 
concentric spheres ; the flux lines are radial straight lines 
diverging from the source. At unit distance from tho 
source the spherical surface has an area of 47r; if tho 
soiirce is of unit strength this is indicated by drawing one 
flux line per unit of area to represent the magnitude of tlic 
vector field at unit distance. A total vector flux of 4W 
lines cuts normally through all spherical surfaces and henco 
the magnitude of the vector field at any point varies with 
the inverse square of tho distance from the source, since 
the spherical area varies directly as the square of the 
distance. If the strength of the source is the flux is 

In Bome fielda the point is an actual source of material flux, 
e.g. when liquid is entering at a point within a fluid ; the vector 
is then tho outward radial velocity of tho fluid. In moat cnaca, 
however, the flux is non-material, e.g. of heat, or of gravitational 
or electric forces. Note that in the case of gravity the point mass 
can only attract another mass placed In its field; the forco is 
essentially radially inward. But In hydromechanics, heat, eloc- 
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tricity, and magnetism both positive and ncgfttlvo flourcca, the latter 
called can exist, making posaiblo both outward and inward 
fluxes. Fluid or heat may cither enter or leave; electric or mag¬ 
netic forces may be cither repulsive or attractive. Massea net in 
tho vectorial senfio in the eamc way as sinks. 

It will be convenient to fix ideas by examining nn actual 
example, namely, the field of electric point charges in 
vacuo. Two such charges of positive electricity hdvc been 
shown by Coulomb’s experiments to be repelled with a 
force hq^q^/r^ in the line joining them. If unit charge is 
defined as repelling a similar charge i cm. distant with a 


S 



Flo. 28.—Electric Field of Point Charge 

force of I dyne, we have the ordinary electrostatic system 
of units, ^ « I and is the force in dynes between 

charges q^ and q^ e.s.u. Let a positive charge q be put 
at O, Fig, 28, then the force per unit positive, charge 
placed nt P is’called the electric force or field strength ; 
its magnitude is q/r^ and its direction that of the unit 
radius vector r^, i.e. of r, then 
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The potential difference between two points A and i? ifl 
defined ns the work done the forces of the field 

in moving the unit positive charge from to B by any 
pnth. Since like charges repel, positive work must be 
done against the field forces to bring the charges nearer 
together. The work tliat would be done by the field forces 
is £‘rfr ; hence the potential difference is, by definition, 



Now let A be at infinity—where 6 is ^ero—and make B. 
coincide with P ; then the integral becomes 



This 13 the poieniial at P, namely, the work done on unit 
positive charge in bringing it from infinity to the point P 
by any path, Using Equation 5,11 8ho\V8 tliat 

£ - ir, = - grnd 

A similar argument can be adapted to other fields. In 
general, let F be a flux vector at a distance r from a source 
of strength q ; then the Newtonian potential is defined 

0 - ! 7 A :.(7.1) 

and F = - grad ^ ^ - grad {q/r\ . . (7-2) 

of which the three cartesian components are respectively 

- - d<f}/dx, ^ d<l>/dy and 7 ^’^ - d<f>/dz . (7.3) 

The vector F points in the direction of the greatest rate 
of decrease of For example, if (f> is temperature it ia 
clear that the flow of heat will be in the direction in which 
the temperature falls most rapidly. 

Since potentials are scalar functions of space they are ' 
arithmetically additive. Thus, if there are sf \ral point 
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sources q^y ... at distances fi, fg . . . from a point 
in the field, ihe potential is 

] • • ■ t?-*) 

Tide may be extended stdl further. If the sources 
are continuously distributed with density a per unit 
area over boundary surfaces the summation becomes the 
integral 



taken over all the surfaces. Each surface element is 
regarded as a point source yielding a Newtonian potential. 

In Equation 5.12 it has been shown that except at a 
point occupied by a source the field F is solenoidal, i.e. 
div F - o ; then 

VV = o:.( 7 - 6 ) 

i.e. Laplace’s equation is satisfied by the potential at all 
source-free points in space, since div grad ^ is, by Equa¬ 
tion 5.12, zero. A field of this kind is, therefore, known 
as a Laplacian field and <f> is called a harmonic function. 
Equations 7.1, 7,4 and 7.5 mny be regarded as solutions 
of Equation 7.6 for various physical conditions. 

3 . Volume Distributions. In many problems the 
field is due to volume distributions of sources, each element 
of which acts according to the inverse square law. Familiar 
examples arc the electric field due to a space charge of 
electrons or of ionized gas molecules, the attraction of 
solid bodies, the flow of internally generated heat (as in 
an electricslly-hcatccl conductor), Sre; Let q be the density 
of sources per unit volume within a closed surface s ; 
(} will be a scalar function of position. Then it dv is a 
volume element in Fig. 29, gdv is an elemental source. 
Two cases arise and will be separately examined. 

6 
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Let P be outside the surface then the total potential 
At P is 



taken over the volume enclosed by s. Also, since qdv is 
constant for each differentiating operation, 

F - - grad <j} and 

exactly as before, The field external to a volume distri¬ 
bution is solenoidal, i.e. the potential is Laplacisn. 



, Fig, 29,—^Potential of a Volume Distribution 
If P is inside the surface s there is a small difficulty, 
since r in the denominator of the integrand can now become 
zero. Surround P by a small spherical surface of radius 
e ; then, taking the volume integral through the apace 
confined between s and the sphere, P is outside this region 
and the integral is finite. Now let the sphere be made 
vanishingly small; then in the limit, its volume appronclics 
which varies with e*. Since r varies witJi dv/f ^ 
depends on and vanishes ns b becomes zero. Ilencc 
Equation 7.7 expresses also the potential inside the suiface , 
if it is regarded as the limiting value reached-when a small ; 
sphere round P is reduced without limit; the field is not, I 
however, Laplacinn, In Equation 6,2 it has been ahowi 1 
that when V<“grad(i/r) the total normal fiux of V 



THE SCALAR POTENTIAL FIELD 


73 


through surface enclosing a source of unit strength is 
« 4jr independently of the position of the source. Hence 
when F - grad 96, where ^ is due to sources of total 
amount no matter how distributed within the surface j, 
the flux is + 471^. In the present case q is the volume 
integral of the source-density, so that 




Applying the divergence theorem, Equation 6.1, 



n'¥ds 



div 



so that div F - 47^^ = - div grad 


and W " .(7'®) 

which is Poisson^s equation. Hence the field in a region 
where there is a volume distidbiitioii of sources is not 
solenoidal. There is divergence) i.e. total normal flux per 
unit volume, at any point, of an amount equal to 4^1 times 
the tlcnsity of sources per unit volume. 

It is obvious that the same results will hold for two or 
more distributed sources, interpreting tlic integral;^Equn- 
tion 7,7, as extending over nil the enclosed volumes. In 
the space between them the field is solenoidal ajul the 
potential Laplacian ; within each volume the field has n 
Poissonian potential and has divergence. 

4 . The Potoiitinl Operation. The integral, Equa¬ 
tion 7,7, taken over nil space may be considered to he the 
solution of Poisson's differential equation, Equation 7,8, 
for any distribution of sources throughout space. 'I’lio 
operation of finding the potential due to a given distribu¬ 
tion of source density q is indicated by the notation pot q, 
due to Giblis. Then, making use of Equations 7.7 and 7.8, 



4JI 


pot VV 


( 7 ''J) 
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Hence the operations 

pot( ) and V*( ) ■ • 

/pi 

Acting on a scalar operand are invdree, one undoing the 
effect of the other in much the same way that integration 
and differentiation are inverse operations. Applications 
of this notation have been much developed by Gibbs, 
and a few simple extensions of the idea will be found in 
the next Chapter. 

5 . Multivalued Potentials, In mathematical text¬ 
books it is shown that the operations gi'ad, div, and cutl, 





Fio. 30.~Multivalued Scalar Potential 


as well as theorems such aa those of Gauss and Stokes, 
are valid provided that the scalar or vector point functions 
to which they are applied remain finite, free from discon¬ 
tinuities and single valued within any given region of 
space and that the same conditions also restrict their 
derivatives in any direction up to the second order (so that 
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operntions such as V ^ have a meaning). Forliinatcly, 
the conditions of finitenesa md continuity are well oh.scrvcd 
in most fields of physics and engineering, but inijxntant 
cases do arise in electromagnetism and hydromcvliariics 
where the idea of scalar potential leads to mnilivalued 
functions. 

Consider first a simple mathematical example, the 
potential given by 

* arc tan {y/cv) “ ± Att, 

^Yhele k is any positive integer including zero, 'fhcfi the 
‘ equipotentiaL surfaces are radial planes 0vnnst.iiit 
iutereecting along the Z axis, shown in Fig, 3 o(< 0 * 
This intersection violates the condition of single-viilucLlness 
mentioned on p. la; the axis is termed a singularity. 
Taking the gradient, 

grad ^ - (^j|. I- j| + k^|)(^arc tan 


“ - 0 I + CC.3 01 ) 

Fhc lines of the vector, gnul nre in the tlircccion of the 
counter-clockwise tangent at imy point P, pcrf>endiciilar 
to the * equipotential * plane through that point. To cal¬ 
culate the line integral of grad (J) we need 





xdy^ ydx 

.V- + y^ 


Consider a circular path round O in the plane of XV; 
then, as a is constant, from ^ 


We liave 


X ^ a cos 0 and y a sin 0 


dx - a sill OdO, dy a cos QdO and xdy - ydx - a-dO ; 
hcnCQ 
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along any portion of a circular path. Since it is inde¬ 
pendent of a the integral applies to cuiy path.* For a 

closed path such aa is clearly zero, since to trace It 

out the angle subtended at O by the area is passer 1 over 
once positively and once negatively by tlie radius vector 
from O* This is true even when A ia infinitesirna! and 
hence curl grad is zero, as required by Equation 5.3. 

Now trace a semicircular path ahy giving a line integral 
of 71 . Likewise, a line integral along the semicircle a'b* 
gives i.e. the same magnitude but of opposite sign. 
But the beginnings and ends of the two paths are the 
same ; hence the value of the line integral depends on 
which path' is followed as well as upon its end points, 
Round any closed path linked about the axis of Z the 
line integral is ± 271 per tour, the upper sign relating to 
rotation in the sense from X to Y. Such a closed path 
cannot be reduced to zero without cutting through the 
Z axis and cannot satisfy Equation 5,3. The potential at 
P, defined as on p. 70, is multi-valued, being increased or 
diminished by m for every tour made round the axis of Z 
before arriving at P from nn infinite distance. 

Mathematicians describe a region of space in which 
these phenomena occur as ' cyclic * about a singularity in 
the region ; in our problem tlie singularity is the axis of 
Z where all the level surfaces intersect. A path such ns A 
which can be contracted to zero without enclosing or 
cutting through the singular point or points is ' re¬ 
ducible ^; a path such as abVa is ' irreducible * since it 
cannot be made zero, but can only be made to fit more 

• Any path can be reaolved into three compononta, viz^ parallel 
to pz, radially from OZ and around it. The first two can add 
nothing to the line integral/eince the linea of flow of gmcl ^ aro 
circular. At any radius a the path around the axis is adO^ for a 
atnall angular rotation. The magnitude of grnd <j> being i/a, 
the element of line integral ia dB, 
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and more closely round tlie aingiilarity. A cyclic region 
can be made ncyclic by inserting impassable barriers to 
prevent the drawing of irreducible paths, thereby isolating 
the singularity. For example, in Fig. 30(a) take the semi- 
iniinite plane of ZX to the right of the Z axis aa a barrier, 
shown shaded. Then 0 cannot now exceed 2 n and ^ is 
thereby made single valued and the conditions of potential 
theory are satisfied. 

Tlie problem just examined ia easily seen to represent tlio 
magnetic field of a current i electromagnetic units in a thin straight 
wire extending to infinity along tho Z axis if wo tako 
^ — 21 arc tun (y/x). 

Then the magnetic force la 

H — — grad tf> ( 2 i 7 n)(*- sin 6 1 -h coa 0 J) 

and the line integral of H round a closed path encircling the current 
ia 4m j tho BO-called magneto motive force. Remove the current, 
niul let the half-plane of ZX to the-right ho occupied by a thin, 
unifonnly-niagnetizod ahcct of moment i per unit urea with tlio 
north-polar surfoco uppermost (Fig. 3o(^)), Then tho mngnetic 
field of tho sheet and the current nro identical In distribution, but 
that of the sheet is acyclic since no path can pass through it to 
encudo the wjrc. Its scokr potential is ^ « 2i(Tr — 0 )^ 0 being 
limited between o and zw (p. 90). 

The fundamental axiom of potential theory is that the 
lino integral of tho vector field is independent of the path 
traversed between two points ; the theory is not nppUcalde, 
therefore, without some artificial aid to problcma \yliich 
do not satisfy this condition. ‘The most imiK>rtnnt ease 
is the essc/Htnlly cyclic magnetic field of nii electric current 
flowing in a linear circuit of any shape. Exploring the 
field, in the usual way, with a unit magnetic polo which is 
moved about from point to point, will result in clianges in 
the amount of flux from the pole which links the circuit, 
inducing electromotive forces therein ; if tho current is to 
remain constant, this will require adjustments to be made 
to the voltage of the battery in the circuit. For every 
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complete tour made by the polo round the current, potcnllal 
energy ± 4 Jit must be credited to the pole; this energy 
is, of course, derived from tlio battery. PJiysically, tlicrc- 
fore, the scalar potential of the field tmai be imiltivalued, 
and its cyclic character has this physical cxpliination, To 
satisfy mathematical requirements, Jio waver, we re [dace 
the field of the current by that of n thin, iiorruaI!y'’njuRnct- 
ized sheet with the circuit as its bouiidaiy eelge-'-the 
equivahut mag7ieiic shell. The potential <Iuo to tJic slicll 
is rco, (o being the solid angle subtended by the circuit at 
any point in the field, and is singlc-vaUied na required by 
potential theory, p. 12 ; the field is then acyclic and the 
problem has been reduced to one in inagiic Costa ties ; see 
further, p. 86. It is, however, not in agreement with the 
essential physical nature of n current field, niul must Iks 
regarded as a mere mathematical convcnicnco, Agrecinciit 
can be secured by treating the field of a linear current fia n 
degenerate case of the vector potential field of n current- 
carrying medium, as on p. 84. 



CHAPTER VIII 


THE VECTOR POTENTIAL FIELD 

*Tho Magnetic Field of n Steady CiUTeiit. 
z^r potential fields occur most commonly in electro- 
letism and in hydromechanics. It will be cojivcnient 
ti oduce the subject by working out the properties of 
ingnetic field of a steady current as a practical ilhistra- 
of tlie fundamental ideas ; the reader will readily 
: them to any analogous physical problem in which 
interested. 

5 have seen on p. 77 that the space round a circuit 
ing a current is a cyclic region ; the physical meaning 
:s 13 that the work done on a unit magnetic pole in 
sorting it once round a constant current of i electco- 
Dtic units is ± ergs. Tliis is the well-known 
theorem, 

laidcr now a medium, such as the interior of a copper 
ctor, which is carrying a current distributed in any 
wny. The current-density at any point is c clcctro- 
itic units per sq, cm, of a surface perpendicular to 
ircction of flow. Draw any closed curve in the 
m and take any open surface having the curve for 
inding edge. If u ds is an clement of vector area on 
rface the current tliroiigli the clement is the scalar 
2t c-n ds, and that through the whole surface is the 
il of this taken over the entire area. Now calculate 
tgential line integral of the magnetic force H round 
sed curve, in the sense of a right-hand screw relative 
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to the unit surface-normal n, as in Fig, 31 ; then the 
circuital theorem makes 



Apply Stokes's theorem, Equation 6.3, to the left side, then 

n-curl H “ 4^*: 

and for this to be true of any surface we must have 

curl H ptc , , . . (8.r) 




open surface s 



Physically this means that we put a very small plane area 
at a point in the medium and turn it about until the line 
integral of H round its boundary is the greatest possible; 
the amount of the integral per unit area is theinagniliido 
of curl H and the positive normal points in the direction 
of the curl. But the equation relates this to the current 
densiCy c, which must, therefore, be passing normally 
through the small area when it is in the orientation giving 
the ^eatest value of boundary line integral. 

It is to be noted that since ^v curl H is zero, Equation 
5.6, div c is zero also. Hence currents are aolcnoidal, 
i.e, flow in closed paths. Further, since there can bo no 
such thing as a space distribution of free magnetic poles, 
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div H is zero also ; the lines of magnetic force must, 
therefore, form closed curves. 

2 . The Vector Potential. On p. 55 we have seen 
that when a vector has curl and no divergence it may be 
represented as the curl of another vector, called the vector 
potential. The reason for this name will now be explained. 
First take at any point in the current-carrying medium, 

H curl A, . . . . , (8.2) 

where A is the vector potential of H. Using Equation 5.7 
on the above, 

curl H curl curl A - grad div A - ^ 4jrc. 

We may, without any loss of generality, assume that A, 
like its related vectors H and c, is eolenoidal; then div A 
is zero and 

^ V^A “ — . (^'3) 

which is a vector form of Poisson's equation relating A 
and c within the medium. Using Equation 5.4, put 
A « + Ayi + AJi and c - + c^k, 

tJien 

VMj - ^ V^Ay = - and 
If t is the distance from an element of volume, where the 
current density is c, to any point in the medium at which 
tlic components of A are to be found, analogy with Equa¬ 
tions 7.8 and 7,9 gives 

(6.4) 

Recombining the components, 

^-potc---potV^A; . (8.5) 

the integrals extend throughout the volume in which c 
is distributed. 

Equations 8.4 show that the components of vector poten- 
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th\ are Crilculated from the components of c in exactly 
the same way as the scalar potential <f> is found from a scalar 
source distribution of density g. Hence the name vector 
potential. Interpreting Equation 8.5, at the i^oint where 
the vector potential is required draw a vector parallel lu 
the current density c at any other point in the mediumj 
the length of the vector being proportional to the ningni- 
tude of c multiplied by dv and divided by r, the distance 
between the points. Repeat for every element of volume 
into which the region containing c can be divided a nit 
sum the vectors obtained ; i.e. integrate their effect at tlia 
point for the whole volume. The resultant vector is the 
vector potential at the point; the same result is more 
prai^cally obtained by calculating Ay, A^ and vectorially 
adding these rectangular components. 

By exact analogy with the scalar potential field, p. ^2, 
the vector potential at a point outside the medium con¬ 
taining c is found by appropriate evaluation of the intcgrnls 
Equations 8.4 and 8.5, but now 

curl H « curl curl A « V s o 

at the point. Each component of A satisfies Laplace's 
equation. H is now a lamellar vector and can be expressed 
as the negative gradient of a multivalued scidar potential 
aay. Hence at points outside a current-carrying medium 
the magnetic force is given by 

H « - grad « curl A. 

At the surface of separation which divides the currcjlt- 
carrying medium from outside space, the two expressions 
for A, in the curl and non-curl regions respectively, reduce 
to equality. Their analytical forma in the two media arc, 
however, quite different.* 

• A VMiety of practical cases will be found in the author'i 
Electromagnetic Problems in Electrical Evgiiieerhig (Oxford 
University Preaa, 1929). 
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3 . The Potential Operation. Tlie quantities A» H, c 
are all solenoldal fluxes and are related one to the other by 
the process of curling, Equations 8,1 and S.2, and the 
calculation of potential, Equation 8,5. Summarizing these 
connexions, 

A - pot c 

H - curl A = curl pot c 

c - — curl H ^ — curl curl A •• —curl curl pot c, 

47r 4^ '4^1 

The fortuitous 1/4?^ is a consequence of the system of 
units that has been used ; it may be abolished by changing 
to the Heaviside rationalized system. The third relation 
shows that the operation denoted by curl curl/4ii: exactly 
annuls tlte operation of computing the vector potential. 
Again, the first and third shows that 

A - pot ^ 

90 tliat pot/4^ and curl curl also annul, Hence it follo^v8 
that pot and curl curl are cofnmutative operations ; we may, 
therefore^ find die vector potential cither before or after 
the operation of double curling. 

Although the three vectors considered are the only 
ones with physical meaning, we can extend the process in 
either direction. Adding a fourth, D ■=- curl C, it follows 
from the construction of the table that H « pot D: Any 
vector in the continuous series is the curl of the one imme¬ 
diately above it (neglecting the 4Jr factor) and the potential 
of the vector two below it. Now 

H - curl pot c - pot D = pot curl c ; . (8.6) 

hence it is immaterial whether wc calculate the vector 
potential before or after calculating the curl. 

These tclations and many others have been developed 
by Gibbs into a complete calculus for manipulating sok 
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noidal vectors, but to pursue it further would be beyond 
our present object. 

4 . Linear Curi'ents, In electrical theory there h 
considerable practical interest in the magnetic field of 
currents flowing in infinitesimally thin conductors forming 
closed circuits of any configuration. Such circuits are 
called Ihiear, since the current is caused to flow along 
straight or curved lines, and are a close approximation to 
practical circuits of thin wire. In hydromechanics, sucli a 
linear curve constitutes n vortex filament in an incom¬ 
pressible fluid, the fluid being in rotation about the lino 
as axis. In both cases the field is cyclic but non-curl. 
The electromagnetic case is often treated magnetostatic- 
ally by the * magnetic shell * conception, but we shall now 
see how vector analysis enables ua to express the field in 
terms of vector potential for a medium free from current 
everywhere, except in the filament constituting the linear 
circuit; precisely similar arguments may be applied to 
vortex filaments in a fluid or to any other linear singularity. 

Let i be the current in the circuit and o' its vanishingly 
small cross-sectional area. If t be unit tangent at any 
point of the circuit, we may put c = it/a and dv a dl^ 
where dl is an element of length of the circuit as shown in 
Fig, 32; then, 

edv ^ itdl i.dl 


and 


A 




will give the vector potential at a point P when the integral 
is evaluated round the closed circuit. The magnetic force 
at P is 


H ^ curl A - j curl - i curl O^dl ^ I'o 


curl ~ dit 
r 


since by Equation 8,6 the operations curl and pot are 
commutative. Using Equation 4.14, substitute i/r and 
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t tJierein; remembering that the curl is calculated at P 
and that t is not a function of that point, ao that curl t 
tliereat is zero, gives 

curl (t/r) - [grad (i/r)] x t - ^ (fi x t)/r^ 
from Equation 5.11, Tj being unit vector drawn from the 
circuit element towards P, Substituting gives 


- to 


U sin 0 dl 


( 8 . 8 ) 

where u is a unit vector in the sense of the vector product 



Fro. 32.—Ampere’s Rule for Magnetic Force of Circuit Element 
(Biot-SgVart I^aw) 

of Fi and t, i.e. t anth u form a right-handed system. 
For example, if and t are in the plane of the paper u 
is perpendicularly outward toward the reader. Equation 
8.8 is known cither as AmpSre*s rule or as the Biot and 
Savarl which states that the magnetic force due to an 
clement dl of a closed circuit has a magnitude isin 6 dl/r^ 
and is normal to the plane containing dl and r ; the direc¬ 
tions of magnetic force, d/and r, forma left-handed system, 
since H is in the sense of — u. 

Calculating curl H, since Equation 5.3 shows curl grad 

* The general rule ia due to Ampfere for any ahape of circuit. 
Riot and Savart^B proof applies only to straight conductora. 
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(i/r) to be zero, the field of linear currents lias no curl. 
As shown on p. 78, scalar potential theory can be applied 
to the cyclic non-curl field by replacing the current by a 
uniformly-magnetized shell having the circuit for its 
bounding edge. The north polar face of the shell is directed 
toward a point P, as in Fig. 33, when the current as viewed 



Fio. 33.—Equivalence of Magnetic Shell and Linear Current 

from P appears to run counter-clockwise. If o) is the solid 
angle subtended by the shell, i.e. by the circuit, at P, then 
tlie scalar potential is taken as 

(f> ^ ico, giving H ^ - grad - tgrad a) . (S.9) 

We now proceed to show that this agrees with the vector 
potential. 

Give the point P a small vector displacement in a 
sense to increase r, the vector joining P to an clcmcJit d\ 
of the circuit; this is equivalent to fixing P and moving the 
circuit bodily in the opposite sense to the dotted positioji. 
The element sweeps out a vector area dl x dp which 
subtends a solid angle Vi ’di x dp/r^ at P, wJicrc Vi is 
unit vector in the sense from P to the element. Tiio 
integral of this taken round the whole circuit is the total 
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change in solid angle subtended by the circuit when it ia 
slightly displaced. Applying the rules for scalar triple 
products, Equation 2.21, gives 

^■(Jl X dp) - X X 

if is unit vector along dp. Dividing by dp and in¬ 
tegrating round the circuit gives 

day fTj' ,, 

since P moves away from the circuit this is the rate of 
decrease of solid angle subtended by the circuit. But 
^rad CO is the total rate of increase of (o ; hence ~ Prgrad (o 
ia the component rate of decrease in any direction denoted 
by Pi, equivalent to dco/dp above. Also if Ti is unit 
vector drawn towards P from ^ 1 , Ti' “ - Ti. Taking the 
scalar product with Pi in Equation 8.9, 

. ,rfa> - 


H'pi - - fpi'grad 

For this to be true of any displacement, 


which is Equation 8.8. Hence the scalar magnetic poten¬ 
tial defined by Equation 8.9 gives the same magnetic field 
aa the vector potential of Equation 8.7; we can write, 
therefore, 

H = curl A - grad 

where A - ^ ^ proving the equivalence of 

the current aiid the shell in the field produced at any point. 

6 . Simple Examples of Vector Potential. The idea 
of vector potential is more readily grasped by considering 
simple examples of linear currents. First take the case 

7 
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of fln infinitely long, thin wire cari’ying current along the 
positive direction of Z as shown in Fig. 34(^2)* The lines 
of inagnetic force are circles with their centres on the Z 
axis, the circles lying in planes parallel to the plane of XY ; 
in other words, since thq wire is infinitely long, the fielcl 
distribution in the plane o£ XY ia precisely similar to tliat 



in any other parallel plane. Since all the elements of 
current are in the same direction, the vector potential at any 
point Tis a vector parallel to OZ and can befound by a simple 
^alar integration, Writing dl^kdz and 
Equation 8.7 gives 


“^1 (*/")] 

-2fk^ arc amh(2^/rt)J 

' + v'(»* + «‘)}} - - 2» log «)k, 
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where la an infinitely great constant arising from the 
infinite extent of the ‘ circuit * round which the integral 
has been taken. Since we are concerned with derivatives 


y 



Fio. 35.—^Vector Potential of Circular Ctirrent 


>f A, the presence of this constant is immaterial; without 
0S8 of generality we can take 

A ^ « - 2fk log rt « - *15 log {x^ + y^), 

3 ince and arc zero, the usual expressioiM for the 
:omponenta of curl A give, from Equation 4.9, 

- 'dAt/dy « - 2 iy/{x^ + - (2i/«) sin 0, 

i/y " - dAJdx - + y^) - (21/rt) cos 

(cf. p 75.) 

Surfaces of constant vector potential are coaxial cylinders ; 
vith equal increments in paasing from one to the next 
heir radii mcreaae in geometric progression. Lines of 
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vector potential nre» therefore^ closely packed together iicnr 
the wire and fall in density at greater distances as siio^vn 
hy 3#)* 

The corresponding multivalued ficalnr potential is cn^ly foUJ» A 
Aa in Fig. 30 (i) take the rlght-hnnd half plniio of ZX^ 
equivalent shell (p. 77). Drnw a epliorc of unit rndiiis with cen ^ 
at P; then the Jialf plane aubtendB an area 4.n(vr — ff)/^7T ^ ^ ^ 

at the spherical surface. This is tlic solid onglc siibtcncic y 
the plane and Equation 8.9 gives ^ 

The components of r- gmd i.e. of H, are identical wK 1 n 
of curl A on p. 89^ ns the reader should verify by using *qiifl- 
tion 4.3 and 0 ^ arc tan (y/x). 

Another simple case is the current in a circle of "'**'*h 
shown in Fig. 35, the circle lying in the pliine of Yifv. 
At any point on the axis of X the vector potential is 
the point being equidistant from all elements of t ic 
circle. For a point off the axis, nearer therefore to one siuo 
of the circle, it is geometrically obvious that the rcsidtiint 
vector potential due to the whole circle is a vector parnJ c 
to the tangent at the nearest part. Let the point describe n 
circle about the axis, parallel to the current ; then by 
symmetry this argument applies to any point on this locus. 
Hence the lines of vector potential for a circular ciurcnt 
filament are parallel circles with their centres on the 
common axis. To give an expression for A and hence fm' 
H requires the use of elliptic functions and is beyond our 
present discussion. 

The diagram, Fig. 35, represents also in one mcridianal 
plane the lines of flow of fluid about a circular vortex 
filament; the complete picture is obtained by revolving 
the diagram about OX, 
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THE ELECTROMAGNETIC FIELD EQUATIONS OF 
MAXWELL 

1 » General. The most striking example of the success 
of vector methods is provided by the theory of the electro¬ 
magnetic held. Although Maxwell introduced the funda¬ 
mental vector ideas into the theory, he did not make any 
use of vector analysis in manipulating his equations, relying 
entirely upon the usual lengthy cartesian methods. It is to 
Heaviside, Lorentz and many Later workers that we must 
turn for the exclusive adoption of the concise vector 
equations and their treatment by purely vectorial methods. 

Our present object is to state the necessary equations 
and to use them as an exercise in the application of prin¬ 
ciples established in earlier cliapters. We shall not, 
therefore, be concerned with questions of a purely physical 
nature, for the discussion of which the reader is referred 
to the special works mentioned in the Bibliography.* 
One main assumption will be made throughout tliis 
chapter, namely, that the medium sustaining the field is 
liomogcncous, i.e. composed of one kind of material in 
the region considered, and also isotropic, i.e. having 
identical physical properties in all directions at every point. 

2 . Maxwell’s Equations. At any instant let c be 
the current density and H the magnetic force at any point 
in the field, both being expressed in electromagnetic units, 

^ Sec particularly F. W. G. White, Electromagnetic Waves 
(Methuen, I934>) in this aeries of Monogrnpha. Also aco p. ii8, 
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Then by the method of p. So the circuital or magneto¬ 
motive force theorem, Equation 8.1 gives 

curl H 4Jrc. 

In deducing this relation, however, a purely cpndiicCing 
medium was assumed, c being the actual flow of electric 
cliarges per second normally through unit arca^ Maxwell 
observed, however, that electromagnetic phenomem take 
place in dielectric non-conductors, even in vacuous space, 
in which electric flow is impossible. To account for tliia 
he postulated a vector D, known as tlic electric dkplacemertt 
at any point and having the physical dimensions of charge 
per unit area, the time rate of change of which is called the 
disphc€in&it current density at the point. He regards tho 
electric force 8, which is the force per unit charge tending 
to move electricity at any point, as related to the dis^ 
placement D in much the same way as stress is related to 
strain in an elastic solid. Using electromagnetic units, if 
K is the dielectric constant of the medium, 


D -/c8/4Jr ..... (9.1) 


In general a semi-conductor will exhibit both phenomenaj 
i.e, true flow of electricity, the conduction current, togctlior 
with displacement or, as an engineer would say, charging 
current. Writing the total current density as 




( 9 »> 


the first equation must be amended to read 


curl H - 47rC = 4jr(c + 6) « 4jrc + wrfi, 

using the dot notation as an abbreviation for time diffei'cn- 
tiation. MaxwelPs assumption that a displacement ciiiTciit 
gives rise to a magnetic force in the same way as a con¬ 
duction current was a stroke of genius that led him to the 
theory of electromagnetic waves in space and the electro¬ 
magnetic nature of light; his conclusions were not veri- 



ELECTROMAGNETIC FIELD EQUATIONS 93 

fied by experiment until many years later. The total 
current flows in a closed circuit, i.e, it is eolenoidal and 

div C - o.(9.3) 


Again draw in the field any closed curve and calculate 
the line integral of the electric force round it. This is the 
electromotive force round the curve and by Stokes’s theorem 
can bo transformed into the surface integral of curl Z taken 
over any open surface with the curve aa its edge. But by 
Farnclay’s law the electromotive force ia the rate at which 
the magnetic flux through the surface is decreasing. Let 
B — he the magnetic induction, or flux density, /x 
being the permeability, then 


0 ^11 n curi 8 di - - ds^ 


which to be true for any surface bounded by the curve 
requires 

curl £ - - 6 - - 

The lines of magnetic induction are necessarily solcnoidal, 
1.0. 

div B - div H - o. 


Summarizing, Maxwell’s equations for the field are 



curl H 4JrC - 4Jic + . 

• • (9-4) 


curl £•^7-6-— 

. . ( 9 - 5 ) 

together with 

C « c + f>, . . . • 

, (9.2 bis) 


div B - 0, 


and 

div D ^ 



Q being the volume density of electric charge at any point. 
These general equations will be returned to in Section 5. 

3 . Conducting Media. When the medium is a con¬ 
ductor there ia no displacement current and the main 
equations become 

curl H “ 4Jtc 


and curl £ - — 
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If y 18 the electrical conductivity of the mediuirij then 
c a y£, by Ohm’s W ; the equations can now be written 
curl B - and curl c - - y 6 . 

Taking the curl of the firsts Equation 5-7 gives 
curl curl B « grad div B - - 4^/1 curl c « — 

But div B ia zero, so that the equation for B is 

.... ( 9 -ft) 

Sirhilnrly, by curling the second equation and noting tliat 
div c vanishes, 

V^c^^jrjuyd ..... (9.7) 

These equations have been solved for many practical cases, 
e.g. eddy current induction in wires, bars or plates. 

They can be conveniently put in terras of vector poten¬ 
tial by defining 

ciirIA«B, . . . . (9.S) 

together with curl B * and curl c *=* - y6. 

Now take the curl of curl A and choose A to make div A 
zero ; then 

V®A « - ..... (9.9) 

The third equation makes 

c-~yA .(9.10) 

Eliminating c 

V^A-4jr/xyA.(O-^O 

4 , Dielectric Media, In a dielectric there can he no 
conduction current, so that Equations 9.4 and 9.5 become 
simply , 

curl H and curl 8 ^ ^ ft. 

Curling the first, 

curl curl H grad div H ^ =- /ccurl t =• 

Since div H o, this can be written 

(v H - o, , , (9.13) 
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with precisely similar forms for B, € and D. The 
bracketed operator is known ns the d^Alembertian, after 
Che French mathematician cVAlembert, who first used a 
particular form of it in the theory of wave motion along 
clastic strings. It is not difficult to show that Equation 
indicates that the components of H (and those of the 
other three quantities also) arc propagated in the medium 
with a finite velocity In free space this is the 

velocity of light, a fact which led Maxwell to the conclusion 
tliat light is itself an electromagnetic phenomenon. 

6. Energy Considerations. Calculation of the energy 
conditions in an electromagnetic field provides an excellent 
example of vector methods. Returning to the general 
case of a semi-conducting dielectric, consider an element of 
volume (b) in which the electric and magnetic forces are 
£ and H respectively and the conduction current density is 
c. Then the Joule heat generated per second within the 
element ia 

dP “ e-c dv s dx dy dz, 

tJie cartesian form enabling the reader to verify the vectorial 
statement. Now Equation 9.4 makes 

— curl H - f), 

and ^ ^ ~ ^ 

Equation 4.12 gives 

div (£ x H) ““ H'curl £ ^ E'curl H, 

ao that JP = r—H curl £ - div (E x H) - E 

|_47r 4^ J 

Using the second main equation, 9.5, and also Equation 9.1, 

dP __i-F/iH'a + X £'t + div (Ex H)] dv- 

4 -’* 
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But it is easy to see that Equation 3.3rt makes 

H-fl - il(H H) - and £• £ - {E«). 

Also Gauss’s divergenco theorem, Equation 6.i, trana- 
forma the volume integral 




the surface integral of the normal component of £ x H 
taken over the surface enclosing a given volume of the 
field, Finally 




H) *1 
( 9 '» 3 ) 


which is Poynt{ng*s theorem. 

On the left side we have the sum of the Joule heat 
generated and the magnetic and electric energies stored 
per second in the volume. On the right is the total 
inward normal flux of a vector 


R - (£ X 

which is known as the PoyfUing radiant vector. It must 
represent the rate of flow of energy per unit area normal 
to R. Note that R is normal to the plane containing 
£ and H and forms a right-handed system with these 
vectors. Equation 9.13 states that the sum of the energy 
per second wasted in heat and that stored electromng- 
netically in any element of volume of an e.m, field lias 
entered that element as a flux of energy per second normally 
through its bounding surface. 



CHAPTER X 


ELEMENTARY PROPERTIES OF THE. LINEAR 
VECTOR FUNCTION 

1 . Tho Linear Vector Function. The simplest 
example of a linear vector function is provided by a vector 
TJ Avlnch is a constant scalar multiple of a second vector V, 
represented by 

.(lo.i) 

k being a constant. This simply means that U is ^ times 
tbo magnitude of V and in the same direction. If V is a 
point function of space, so also is U but altered in scale k 
times. Familiar examples are the relations between B 
and H or between D and E in homogeneous, isotropic 
media. 

Rather more generally, let each component of V be 
multiplied linearly by a constant, different for each com¬ 
ponent. Using rectangular co-ordinates, if 

V-F^l+iy + nk . . . (lo.z) 

And U^i + U^i + UJc, . . . (10.3) 

then this more general type of linear vector function 
requires 

“ KVx. Uy " h,Vy and U, - KV, . (10.4) 
The vector U is not only different in magnitude but also 
different in direction from the vector V. As e^mples 
may be mentioned the relation behyeen D and E in cryat^s, 
or that between stress and strain in a crystalline solid ; 
indeed any case of anisotropic media In which the rect- 
'angular axes htg prwdpal axes of strain. 

IN'Tost general of all is the case of complete anisotropy, 
97 
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where eoch component of U is a linear function of all three 
components of V. In rectangular co-ordinates the relation 
between the components can be written 

The nine Coefficients characterize the transformation ol 
the components of one vector into those of the other ; m 
general, hap transforms Vp into one of the three parts of 
^ay giving a rule whereby to remember the signifienneo of 
the subscripts. The final result of the transformation is 

U - (PV,.(10.6) 

which may be regarded as a generalization of Equation 
10.1 ; 0 is an operator turning V into U by a linear modi¬ 
fication of components specified by Equation 10.5, It ia 
rather graphically described as a Cartesian ie?isory or simply 
as a ie/tso 7 ‘t since the relation of strain to stress in an aniso¬ 
tropic elastic solid obeying Hookers Jaw is of this kind. 

The operator ^ ia also sometimes called an ajfimr, for tho 
following reason. Let V have tho particular value of i*, the 
radius vector from the origin to n point in space. Then U is n 
vector to a second point related to the first by a general linear 
strain ; the points arc said to be in i : i corrcspondenco or in 
affinity. 

2 . Simple Types of Tensors. The essential part of 
fl tensor opeif^tion ia the array of coefficients, such as 
which in the general case are nine in number ; this can be 
symbolically represented, writing tho array in the form of 
a matrix thus 



• More correctly a square matrix. The general type of matrix 
has mu elements in m rows of n columns. It rep resents the nrrny 
of coefficients in m linear equations with n variables and ia of a 
much more general typo than we have apace to examine here. 
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Tenaorfl and their matricea are manipulated by certain 
rules, ^vhich are chosen to give results consistent with the 
algebra of groups of linear equations, such as Equation 10.5. 
A matrix is a convenient symbol for a whole group and 
enables the set of equations to be handled as a single 
entity. It is beyond our purpose to deal in any detail with 
matrix algebra, but one or two simple rules will be 
stated. First, if V « A + B we define 
45 (A + B) - + CPB ; 

i.e, the tensor operation is regarded as following the same 
of distribution as for multiplication in scalar algebra. 
Second, if and 0 ^ are two different tensors, 

^iV + + 0 ,Y = ( 0 , + 02 )y ; 

i.e, tensors follow the commutative law for addition, Tliia 
tlie reader may readily verify by setting down the sum of 
t\70 sets of equations like Equation 10.5, when it will be 
seen tliat the sum of two matrices is a new matrix in which 
the elements are the sum of the elements of 0 i and 0 ^. 
Third, a tensor is regarded as negative when every coeffi¬ 
cient in its matrix is reversed in sign. 

If the columns and rows of a matrix are interchanged, 
the resulting tensor is the conjugate tensorf for which the 
matrix is 

rKx kffj, 

When two tensors are such that one is the conjugate of 
the other, inspection of Equations 10.7 and 10.8 shows that 

Kv " " ^' 1 ' • ■ 

If a single tensor satisfies this condition it is called 537/1- 
TNeirical or self-conjugate. It has only r/A’ independent 
elements and may be written 

Our work is limited to the cartesian tensor and ita square matrix, 
for which w — — 3, 
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“ hy Ku Kt • • • 

^ys 

Symmetrical tenaors arc of frequent occurrence in physical 
problems and their meaning will be adverted to in Section 3. 

In particular when the relation between U and V is 
referred to principal axes all coeflicients arc zero except 
the diagonal corresponding to Equation 104. 

Wlien two tensors are such that one is minus the con¬ 
jugate of the other, the coefficients must be such that 

^nd ky^^^'—kiy (10,11) 
If a single tensor satisfies this condition it is called anti’' 
symmeiricaly anii-self-conjugaie or shew^ It has only three 
coefficients and its matrix is 

po k^r\ 

/{^, . . .(10.12) 

L- lt„ - k„, o J 

Since a skew tensor shares with a vector the property of 
having only three components, the operation of (Pg* on n 
vector V ia exactly equivalent to the vector product of 
two vectors; since the final result is itself a vector U. 

3. The Symmetrical Tensor. The symmetrical 
tensor, Equation 10.10, has a simple geometric interpreta¬ 
tion that will now be explained. Let i, J, k in Fig. 36 be 
unit vectors along rectangular axes Xy Y, Z. With the 
same origin let a, b, c be unit vectors along a second set 
of rectangular axes i 4 , 5, G. The two systems arc related 
by the cosines of the nine angles between all possible 
pairs of unit vectors, one in each system ; these nine 
cosines may be represented by the matrix, 

i J k 
h /jj-n 

b Wli Vit W>a 

cLmi J 
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e,g. 18 the cosine of the angle between i and a, &c. Since 
both systems arc rectangular, it is clear that only tliree of 
the cosines can be independent, i.e, one system can be put 
into the position of the other by tlirce independent angular 
movements* Project the 1 , J, k system into the directions 
of a, b, c; then it is geometrically obvious that 

a - /li + /J + /,k 

b + WjJ + Wfik 

c “ + ?<ok. 


K 



Now from Equation 2.4, = c® ^ i and from 

Equation 2,2 a-b ^ b-c ^ c^a « o ; also ^ k“ - i 

and l‘J - j 'k - k-i « o. Calculating these scalar products 
gives 

+ /a® + /a** - + «8* "■ 

Ijfn^ + /gWia + ^ ” 

Hill + Wa^a 

which are the well-known six relations between the cosines ; 
showing that only three can be independently chosen. 

If the relation of U to V is given in one set of axes, 





102 


VECTOR ANALYSIS 


what is it with reference to the other ? Let the a, b, C 
ayatem be selected in such a way that they arc principal 
axes ; thcn» as defined by Equation 10.4, 

C/a “ KV,, - hVt, and U, = k,V, 
is the relation between the components along B and C 
respectively. Now resolve Vy^ F,, the components of 
V in the X, Y, Z system, into the other system, thus 

Va " /iF, + hVy + hV, ^ t/, ^ kaVa 

Fft “ mjF^ + m^Vy + 7W3F2 \ and - ^/.Fj 

F, - «,F^ + n,Vy + J 

Now resolve 17 ^, C/,, upon A', giving 

- /iC/« + « (Vi" + 

Likewise resolving upon Y and Z, 

* (KhU + 

+ [kJi^ + 

+ 4 * k^1hns)Vt) 

- (Vi^ + h^niin^ + kfpiin^)Vx 

+(Aa^a^ 34 *ft 67 «iW 2 a+fto 7 ;a 7 ; 8 )Fy+(Aa/ 3 H/e 67 / 7 a^+M 8 ^)I^** 
Comparing these three relations with Equations 10.g it 
will be found that the coefficients satisfy Equation 10.9 ; 
the tensor relating U tp V in the system of axes X, Y, Z 
is, therefore, symmetrical like Equation 10.10. It follo^VB 
that any symmetrical tensor corresponds with a transfor¬ 
mation froan tlie principal axes to another rectangLilair 
system of reference. 

It is possible to give a simple graphical construction for 
a symmetrical linear function. Using Equation 10.10 the 
function is 

U - ^ V 

Now project V upon U, i.e. calculate by Equation a. 10 
VU - + VyUy + V,U,. 
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Using the ayminetry of Equation 10.10 in Equation 10.5, 
V-U - + A„F.a + 

+ ky,V„V, + a S, 

aay, where iS is a scalar point-function. But ' 

Ux “ kxxVx + ItxyVy + ]txiV, “ ^ gp- 
and similarly C7„ - J ^ and Uj ^ 

oyjf oyg 

But theae are the componenta of J grad S' in terms of 
Fjj., Vy and as co-ordinates, Equation 4,3, so that 

U « + C/^j .h C/,k - I grad S ; . (10.13) 



Fio. 37.—^Thc Tensor Ellipsoid 


it follows, therefore, that U is a vector perpendicular to 
the surface S ^ cofisiaiii in the direction of the outward 
normal. But *5 ^ co 7 istant is an equation of the second 
degree in the rectangular components of V ; regarding 
these as co-oKlinates defining the extremity P of the 
vector V, the locus of that extremity is a quadric surface 
or conicoid. In particular, if *5 - i the surface is, under 
certain conditions, the teusor ellipsoid, and 

V^U ^ (resolute of V in direction of U)f/ - i (10.14) 

Hence this resolute in the direction of grad S is the recip¬ 
rocal of the magnitude V of the vector U. Fig, 37 shows 
8 
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this geometric relation for a tensor ellipsoid set in the direc¬ 
tion of its principal axes. 

4 , Resolution of a Tensor. Return no^v to the 
general tensor, Equation 10.7, and apply the distributive 
law to it and its conjugate ; then 

+ + . . (lo.is) 

If the reader will write down the matrices for cP + 0 ^ and 
0 - 0^ by adding corresponding coefficients, it is easy to 
see that the former is a symmetrical and the latter a skew 
tensor. Hence any tensor relation between two vectors 
is equivalent to the sum of a symmetrical relation and n 
skew relation. The first corresponds to a change of axes 
without deformation of the system of reference, as has been 
shown, and the second to a. transformation of the angles 
between the axes themselves, 

5 . Repeated Tensor Operations. In vector tlieory 
we frequently need to express relations between more than 
two linearly connected vectors in an anisotropic medium. 
Thus in electromagnetism there is a tensor relation between 
H and B, and a second relation between B and the ponclero- 
motive force on magnetized solids. In a general way, if 


U - 0 iV where 0 ^ 

and V " 0 aW where 0 , 
what is the meaning of 


-B 

B 



U - 010,W } 


The answer is easily found by writing out the thi'eo rela¬ 
tions between the components of U and those of V, Equa¬ 
tion 10.5, and the analogous set giving the connexion 
between the components of V and those of W. From the 
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second set of equations substitutii the components of V 
in tile first set, giving U in terms of W. The final matrix 
for the complete transformation can. be verified by the 
reader to be 

r ^ (^xai^xy+(^a‘xCzi+ ^icyCyt 

1 \I 

L(«ftff’rj+ hzyCyx + ffzz<^zx) + ^zy^yy ++^X2^rz)_) 

. (10.16) 

An inspection of this apparently complex matrix will show 
that its construction follows u very simple rule : The 
element in the rth horizontal row and cth. vertical column 
of is obtained by multiplying each element in the 

rth row of by the corresponding element in the. cth 
column of 0 ^ and adding the results. It is iisual for¬ 
mally to define 010^ as the * product * of the two matrices, 
the rule for the .* multiplication ' being that just stated. 

It 13 obvious that 0 ^ 0 i will, in general, be quite different 
from 0 i 0 ^\ hence the rule for matrix ' multiplication *, 
as also for successive tensor operations, is not commutative. 
The algebra of tensors and their matrices, therefore, fails 
in the commutative law for multiplication exactly as wo 
liave eiirlicr seen to be true of vectors. 

A matrix of particular interest is 

Fi o 6^ 
u- o I o , 

L° ® 

which, as the reader can verify, has the important property 
that 

a 01 ) ^ 0 , 

Hence v has the same effect in tensor algebra as unity in 
ordinary scalar algebra, since it leaves a vector unchanged, 
and is called the imit tetisor or imit 7Jiairix. 

64 The Dyadic. Although tensors are widely used to 

8 * 
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express linear vector relationships in an anisotropic fields 
there is another method by 'which this can be done. The 
dyadic, introduced by Gibbs as an essential part of his 
treatment of vector analysis, is preferred by some writers, 
and its use will now be briefly explained. 

Returning to Equation 10.3 and inserting the values of 
Uy and given by Equation 10.5, the expression for 
U can be written as 

say. Using Equation 2.9, - i*V, Kj, - J'V and 7 , « k-V, 

giving 

Uc.U-V+mj'V+iik-V-(ll+mJ+nk)'V-lP'V . . (10.17) 

The bracketed expression is known as a dyadic and consists 
of the sum of three terms, li, mj, nk, which are its dyads. 
Each dyad is made up of two vectors in Buccession, the 
first vector being called the aniecedmit of the dyad (such 
as 1, m, n) and the second vector its comequent (such as 
i, J, k). It is to be noted that a dyad is not a simple product^ 
but merely a symbolic juxtaposition of two vectors. In 
interpreting the operator the vector nearest to V in any 
dyad is involved in a scalar product with V, this product 
then being a scalar multiplier for the remaining vector of 
the dyad. Comparing Equation 10.17 
that the operation of a dyadic on a vector by means of this 
scalar product rule is equivalent to the^ operation with n 
cartesian tensor. 

In Equation 10.17, 've interchange ]F and V, 

U'-V li+V mj 4 -V.uk=V (li+mj+nk)-VV^. . (10.r8) 
which is a linear vector function different from U and called 
its conjugate. In Equation 10.17 dyadic appears as a 
prefactor and in Equation 10.18 as a postfactor^ 

Starting again from 

U - Ea-l + P^j,m 4 - 
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find commutating the scalar products for Vai Vyt gives 
the identity 

U-Vll+VJm+V-kn-V(il+Jm+kn)-V*V^o^Y^‘V( 10.19) 

The dyadic 1/^^, in which the r 61 cs of antecedent and con¬ 
sequent are interchanged, is the conjugate dyadic of IP, 
Hence any dyadic used as a prefactor is equivalent to its 
conjugate used aa a postfactor. If a dyadic is such that 
W y - it is self-conjugate ; if - - V IP it is 
antiself-conjugate. These definitions are exactly parallel 
with those for the corresponding tensors on p, 100. 

A dyadic can easily be expressed in terms of the nine 
possible dyads of the unit vectors i, j, k, and in this so-called 
tiouio}i form the relation bet\veen dyadic, tensor and 
matrix are even more clearly displayed. On p. 106 substi¬ 
tute in U for Vy, in terms of the unit vectors, then 

+ (VJ + * Ji + 

+ {k^ik + + ft„kk) , . (10.20) 

and also 

+(ft„jkl+/{„kj+/:„kk)] + (/[,^k+ftjJk+A„kk}] (10.21) 

wliifh arc exactly analogous to the matrices, Equations 
10,7 and 10,8. 

If a dyadic is such that when applied to a vector, either as 
a prefactor or as a postfactor, it gives the vector itself, it is 
called the imit dyadic or idemfactor* For this to be the 
case all coefficients for which a +must vanish and 
^xx “ ^ Kz ^ I ‘ Then, ns is obvious, 

J - 11 + Jj + kk makes J-V - V*/ s V, . . (10.22) 

which is equivalent to the unit matrix, since this also 
corresponds with U ^ V. The idemfactor and unit 
matrix occupy the same place in dyadic algebra and tensor 
calculus as unity in ordinary scalar calculations. 
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It is outside the scope of this book to enter into further 
details of the dyadic algebra developed by Gibbs or that 
of matrices used in dealing with tensors. It will bo 
sufficient to say that both form complete, parallel and 
equivalent systems, and have been applied to a wide 
variety of problems in pure and applied mathemniics. 
They are of the non-commutative class in regard to certain 
products which, as with vectors, are matters of dcfiniti(jn. 

7. Application of Linear Vector Functions. It is 
a matter for personal preference whether one uses the 
tensor or dyadic notation in practice, though present 
tendencies would appear to show that tensors and matrices 
are nowadays the more usuah The simple cartesinn 
tensors we have considered are fundamental in all prob¬ 
lems of anisotropy, whether magnetic, electric, thermal or 
elastic, and also in certain purely geometric applications to 
kinematics and kinetics where transformations of rect¬ 
angular axes are involved,'as for example in the treatment 
of finite rotations. • 

More complex tensors, which we have no spaco t6 
discuss, occur in the transformation of spncc-timc co¬ 
ordinates from one obsei'ver to another in the theory of 
relativity. Also, in the general theory of electromagnetic 
machinery the properties of a number of electiic circutlB 
in relative motion are connected by sets of linear equations: 
The use of tensor and matrix notation is a natural means 
for investigating such problems as G.ibriel Kron baa 
recently shown with great success, but this very special 
application is entirely outside the purely vectorial scope 
of this volume and belongs properly to the general algebra 
of tensors. 
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Tins most important operations of vector analysis have been 
expressed in the text in terms of components referred to rect¬ 
angular or cartesian co-ordinntes* Summarizing, they are 

.. ( 4 - 3 ) 


grad VS-—i . 

OX OJ/ dff 


divV=VV = ^* + ^'+^-. 






OAl" 


. ( 4 - 6 , 4 - 8 ) 

(4,9, 4.10) 

. . (5.1) 


iS is n scalar and V a vector point function. These usually BiifTicc, 
but I a some problems it is often more convenient to use soma 
other system of co-ordinates specially adapted to the matter in 
band ; two of these will now be examined. 

When there ib symmetry about an axis. Such as 02 , the 
cylindrical polar or columnar co-ordinates showi^ in Fig. 38 
arc appropriate. The position of a point, P is specified by the 
polar co-ordinates r, 0 and the axial distance these being 
related to the cartesian system by 

X cos 0 ^ y ^ r sin 0 and ar ■» «. 

It is possible to transform the above operations from ono syatem 
to the other by purely mathematical mej\p3, but it is more con¬ 
vincing to work out the operations afresh for the element of 
volume dr, rdO- d;0f shown to a larger scale in the lower part of 
the diagram, (r, 5 , ») being the co-ordinates of ita mid-point. 
As an example, if Vr be the radial component of a vector V at 
F the flux through the element in the r direction is 

( F, + 1^'*) {»■ + yr)<iO da - ( r^r ~ 

which to the third order of small quantities is 

dO dr d» + V,d 9 dr da. 
dr 
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Since the volume of tho element is t dO dr dz, the rndini com¬ 
ponent of div V is 


^r r 




The other components are found by a similar process. Like¬ 
wise, by the method of pp. 43-5 the componente of curl V can 


y 



bo calculated ; , this Is left to the reader ofi* nn oxcrciso. Let 

a, b, k be unit vectors nt P, a being in the,positive direction 
of r, b in the sense of the positive tangent to the circle of radius 
r at P (l.e. along the arc rdd) and k along OZ. Noto that n, 

b, k form a mutually perpendicular, right-handed syatem ; n 
and b, though of unit magnitude, are functions of the poBition 
of P. Then, if Vrt Vgt Vg are components of tbe vector point 
function V in the a, b and h directions respectively, it is easy 
to show that 
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grad 


or r bo 

0* 




divV = 

V V 


I bPi) 
r bfl + 

hVn 
bjS * 



curl 

V- 

V X 

\r bff 

bFA 

f)a/“ 


2>r , 





1 / i)S\ 


•KS 

(rVe) ■ 

i>V,\ 
bO ) 

cliv grad 



, 1 

r* &e» 

b»«- 




K 



When there is fiymmch7 about a point, as in many geometrical 
and dynamical problems, the spherical polar co-ordinatea sho-wn 
in Fig, 39 are useful. The position of P is given by the radius 
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r, the angle of co*Iatitudc bchveien OZ and r, and the angle of 
longitude ^ betAvecn the ZOX plane and that which contains 
OZ nnd r. The relation to the cartesian system is 

a ; r sin Q cos y ^ t sin Q sin ^ and jot *= r cos 
Taking0, <ff) as the co-ordinates of the middle of tlic volume 
element r sin let n, b, c be a right-handed orthogonal 

set of unit vectors at P in the positive directions of these in ere- 
hienta ; note that a, b, c are functions of the position of P. 
If Vrt Vsy are components of o vector point function V in 
these directions, it' is easy to prove for. the volume clement con¬ 
sidered that • 


grad iS « V 5 


1 DP I DP 

■?Jf 7-pineal ■ 
t a. 


r sin 0 dJ * 


di, V - v-v - t 

™,I V . V X y - .in«) - 'g']. 

di.p.dS . V'S 0^) 


I 



PROPERTIES OF V AS A FORMAL VECTOR 


On p. 33 a definition of V in cartesian vector fomi has been 
given, but it must be remembered that V cannot be regarded 
as Q true vector since b/dy, ^/bjar are differential operators 

and not true scalars. The symbol V actually denotes an in¬ 
variant vector differentiator and, as auch, requires an operand, 
■whether scalar or vector, upon which it can appropriately act. 
Its expression in vectorial form, therefore, can only be taken as 
iymholic or /omfl/, inasmuch as V behaves like a vector in formal 
asBocIative and distributive differentiations such aa 


V(t/+ IF)- Vtf+ VIF 
V(17R0-(Vi:/)TF+ 


>vhere XJ and W nre scalar point functions ; these and many other 
simple openUiona the reader will find in die text-books and may 
easily verify. 

It must not be assumed, however, that V In more complex 
operations nccessoyily obeys all the rules of vector algebra, since 
V is only fonnally and not actually a vector. Caro must be 
taken to see that rigid application of those rules does not d«troy 
tile differentiative effect of the complete operation. Particular 
care is needed in using the commutative law; examples arc 
found in the operations V’V and V xV on a vector point 
function V, which become meaningless when commutated (see 
pp. 42 and 45). ’ Wo conclude, therefore, that V ia essentially 
and primarily an operator ; secondarily, and then only os an 
analytical convenience, is it expressed in cartesian vector form. 
In arranging more complex operations the object should always 
be to choose auch a use of the vector laws as will leave differentia- 
tive operations available to act on the operand. 

An important practical case is that of the operator defined by 


C'V where 

C « C«,l + CJ + C^k 

is a constant vector, i.e. one independent of x, y, ar. Taking the 
Bcalar product gives 
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which is 0 Bcalmr different inter applicable either to a scnlnr or lo a 
vector point function. Note thnt V-C would vanish and thus 
have no meaning. In particular let 

c « + c J + Cak 

be a vector of unit length in a given fixed direction defined by 
^ibj its components along the axes, i.c. by tlio diiecllon- 

coaines of c. Then the operation 


^ h h d 

Apply this operntor to a acalnr point function then 


(c-V )5 


55 55 


Now from Equation 4,3 


. ^ „„ 55 b5 

grad 5- V5-:^l + —j 
ox 531*' 



IB the greatest rate of increase of 5 at o point in mognltudo and 
direction; and C'(grad 5 ) is the resolved or component mngni- 
^de of the rate of increase in the direction of c, a scalar point 
function. Calculating this scalar product and comparing results 
shows that 

(c.v)5- c-(grad 5) = c (V5). 

This is known as.the direciioftol derivaiwe of 5 in tho direction 
of c, exprewing the magnitude of the rate of increase of S in this 
specified direction at any point. Similarly, for n vector point 
function 


V . 
(c*v)V . 

where (C‘V)l^a = 


(c*v)E^ + 
dVa dV^ 
''■‘dx 






is Ae diTj^tional derivative of V in the direction of cit ia r 
vector pomt function. 

This theorem provides an alternative proof of a result giver 
on p. 47 , namely, to calculate V x o> x r, where 

to “ <ti„l 4- ttij 4 - 

is a constant vector and 


r - id + 3> j 4 - ;Qrk 
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a variable one. Expanding as n vector triple product by 
quntion 2,22^ 

V X o> X r = (V r)co - (V* 03 )r 

Lit 


/ 5 5 ^ \ . . , .. bV 53 

“ (fe + ^ + s " 3- 

ow from above V*o> baa no opertotional meaning ; we Interpret 
, therefore, ns C0‘V. Substituting to for C, 


id 


5 5 5 

tO'V ~ o,^ + <0^ + O,^. 

(w V)r - ("I^. + yj + 

■" <11^ + + <Oik ^ (o. 

inally, 

Vxcoxr -300-0)“ 20), 
i found by direct expansion on p, 47^ 
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obtain many useful ideas from Vol, 1 of EUcironiagnetic Theory, 
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NOTATION 

■iIB notation uaed in this book is thnt of Gibbs, with the additional 
ioption from Lorentz of tho tenn ^ grad * to represent the 
LITercntiatton of n scalar by the operator V - On the Continent 
setain differences of notation occur and these niso appear in some 
ngllsh'books (often translations). These are (i) the uae of 
.B instead of A-B, or the employment of round brackels (AB) 

) denote the scalar product, (ii) The use of [AB] for the vector 
roduct AxB. (iii) The substitution of ‘ rot ' for ‘ cud *. (iv) 
"howritingof A.Borof A; B for the dyad AB. These difTercnccs 
re easily remembered and vector equations can be read with 
qual facility eitlier in tho Gibbs or the Continental system after 

Uttlc practice. . 

Heaviside uses an unaymmetricol notation that makes his 
writings unnecesssrily hard to follow \ ho usea AB and VAU for 
he scalar and vector products respectively. 

In reading vector equations, expressions such as A B and 
^ X B aie read ‘A dot B* and * A cross B* respectively, 
jiinilarly, is read ^ del S *, V*V is * del dot V ' and V x V 
a ‘ del cross V 

The notation for line, surface and volume integral^ follows 

Iribbs, i.c. 

J( )dl, JJ( )* ,nd jn( 

■cflpcctively, the number of integral signs being a useful reminder 
:hat dl has one, ds two and dv three dimensions of length srresiicc- 
tivQ of any particular co-ordinate system. 

Historically tlie terms used in vector analysis arc of great 
interest. Although the operator V "was introduced by ll.imiJioit 
ho left its development to his disciple P. G. Tnit, who adopted 
the name ' nabla * for this symbol, following a suggestion of 
W, Robertson Smith. Mftx^vcll suggested ' ailed but the 
name * del ' due to Gibbs is preferred for shortncii an<l euphony. 
To Maxwell are due the terms * slope ’ (for * grad * conver- 
gciico ' (for - ‘ div ’), ‘ curl ’ and ‘ concentration * (for V ‘j. 



NOTE ON MAXWELL’S EQUATIONS 


Thu deduction of Maxwell’s equations for tlie electro- 
magnetic field given on p. 92 ct seq. follows thnt given in 
Chap. IX, Vol. 11 of hia Treatise on Electficiiy and 
Ma^neiim, all quantities being expressed in electromagnetic 
unite. It is very common practice, however, to state c, 
H, fj, and B in the electromagnetic system, but £, k and D 
in the electrostatic system. If v is the velocity of light 
m free space the last three quantities are expressed iu 
electromagnetic units by writing v C, and D/v* 
Equation 9.2 then becomes 


V J^nv 

and the field Equations 9.4 and 9.5 are 


curl H - - 4jrc + - £ 

V 


curl C 


together with 


I 

V 


6 - - 

V * 


div B - O 


and div D - p 

The reader will find it a useful exercise to work out nnv 
modtficatbns ntroduCed infn T7„ ^ ^ 

the use of thic • . Equations g.6 to 0,13 by 

ttie use of this mixed system of units J ^ 
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ddition of vcctora, 6 
ffinor^ 98 

A 1 embertian operator, 95 
mpftrc’a rule, 85 
ngular velocity, 46 
compounding of, 46 
relation to curl, 47^ nj 
saocifltive law, for ecalare, 3 
for vector sum, 8 

iot-Savart law, 85 

Sircuital theorem, 79 
lommutative law, for acalara, 3 
for scalar product of two 
vectors, 15 
'url of curl, 53 
lurl of gradient, 51 
iurl of vector field, 42 
in rectangular co-ordinates, 45 
na operation V x> 45 
examples of, 45 

rotation of rigid body, 46, 

motion of fluid, 47 
electromagnetism, 48 
application of Siokea’a 
theorem to, 64, 65 
invariance of, 66 
in cylindrical co-ordinatca, 
III 

in spherical co-ordinates, 11:^ 
Cyclic regions, 76 
Cylindrical polar co-ordinates, 

109 


Decomposition of field into 
lamellar and aolenoidal 
components, 55 
Del operator, 33, 113 
Differentiation of vectois, 29 
sums and products, 30 
triple products, 31 
partial, 32 

del operator, 33, 113 
Distributive law, for scalara, 3 
for scalar product, 17 
for vector product, 24 
Divergence of curl, 53 
Divergence of gradient, 51 
Divergence of vector field, 39 
in rectangular co-ordinates, 40 
as operation, V*j 4a 
meaning of, 41 
Gauss's theorem, 59 
invariance of, 65 
incylindricnl co-ordinates,! x i 
in spherical co-ordinates, 112 
Dyadics, 105 
dyads of, 106 
conjugate of, to6 
nonion form of, 107 
relation to tensor and matrix, 
107 

Fields, scalar, i2 

level surfaces in, I 3 
vector, iHr 

lines of flow in, 13 
Finite rotations, 7, 46 
Flux, line of, I2 
total, 20 
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GauBs’e theorem, 59 

applied to inverac aqimro law, 
61 

Gradient of divergcnco, 52 
Gradient of scalar field, 33, 35 
ns operator V, 35 
meaning of, 35, 38 
in rectangular co-ordinntca, 
36 

in cylindrical co-ordinates,! 11 
in aphericfll co-ordinates, 112 

Helmholtz*a theorem, 55 

Idem factor, 107 
Inverse square law, 57, 68 
point sources, 56 
Gausses theorem, 61 
applied to electric field, 69 
Irrotational motion, 48 

Lamellar vector, 37, 38, 54 
closed line integral of^, 38 
Laplace^s equation, 54, 71 
Laplace's operator, 51 
as div grad, 51 
in rectangular co-ordinates, 
SI 

1 n cy 1 ind ricnl co-ord ina tea, 111 
in spherical co-ordinates, 112 
Level surface, 12, 34, 67 
Line integral of vector, 18 
meaning of, 19 
in lamellar field, 37. 38 • 
Line of flow, 12 
Line, vector equation of, 7 
Linear currents, theory' of, 84 
equivalent magnetic shell, 86 
Linear vector function, 97 
principal axes, 97 
cartesian tensor of, 98 
conjugate of, 99, 106 
skew, 100 
symmetrical, 100 


Linear vector function —cotttd 
resolution of, 104 
repeated, 104 
dyadic of, 107 . 

nonion form of, 107 

Magnetic field, of straight wire* 
77 

of current-carrying tncdhim, 
79 

of linear currents, 84 
Magnetic shell, 78 

equivalent to linear current 
86, 90 
Matrix, 98 
conjugate, 99 

skew and symmetrical, 100 
resolution of, T04 
product, 105 
unit, 105 

Maxwell's equations foi* e.m, 
field, 91, 118 
in conductors, 93 
in dielectrics, 94 
energy relations, 95 
Multivalued potentials, 75, 82 
cyclic and acyclic regions, 76 

Nabla operator, 33 
Newtonian potential, 58 
Normal, positive sign of, iQ 
Notation, comparison of, 117 

Operator, del, 33, 113 
35 . 36 
div, 40, 42 
curl, 45 

div grad, 51, 52 
grad div, 52 
div curl, 53 
curl curl, 53 
potentml, 73, 83 
d'Alcmbertian, 95 
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Per tin 1 differentiation of vcc« 
tors, 3a 

Plono^ vector equation of, 17 
Point function* iz 
Point Bourcc, 56, 68 
PoiBson’fi equation, 54, 73 
Polar co-oi'dinntcs, 109 
Potential operator, 73 
Potential, scalar, 38, 67 
of point source, 57, 68 
Newtonian, 58 
in oicctric field, 70 
of surface distribution, 71 
of volume distributions, 7a 
multivalued* 74 
Potential, vector* 55, 81 
in lamellar region, 82 
of linear currents, 84, 90 
Poynting’s theorem, 96 
Products more than three 
vectors* 27 

Productfl of three vectors* 25 
Products of two vectors, scalar, 

vector, 20 
PseudoBcalara, 26 

Rotation* 43, 48 

Scalar algebra, 2 
Scalar, definition of, i 
Scalar potential fields* la, 38 
of point source, 57* 6S 
level surfaces in* 67 
Inmellnr nature of, 68 
of surface distributions, 71 
of volume distributions, 71 
multivalued, 74 

Scalar product of two vectors, 

14 

commutative property of, 15 
conditions for zero* 15 
Relf, IS 

application to unit vectors, 16 


Scalar product of two vectora 
—cofitcL 

distributive law, 16 

inrcctangularco-ordinatesp t 7 
differentiation of* 31 
Scalar triple product of vectora, 

in rectangular co-ordinates, 
26 

Solcnoidal vector, 41, 54 
Solid angle* 61 
and magnetic shell, 78 
Spherical polar co-ordinates,! 11 
Stokca^s theorem, 63 
Subtraction of vectors* 6 ' 
Surface integral of vector, 19 
meaning of, 20 

Tensors, cartesian, 98 
conjugate, 99 
skew, 100 
symmetrical, 101 
representation by ellipsoid, 
103 

resolution of* 104 
product of, 104 
unit* 105 

Triple product, scalar, 25 
vector* 27 

differentiation of, 31 

Unit dyadic, 107 
Unit tensor, 105 
Unit vectors* 10, i6, 22 

Vector area, 22 

of closed surface, 23 
Vector, definition of, 2 

graphical represent0tion of, 4 
scalar multiplication, 5 
polar, 5 
axial, 6 

Vector fields, 12 
classification of, 54 
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Vector potential, 55, 79, 81 
calculation of^ 8a 
in lamellar region, 8a 
of linear currentB, 84 
of straight wire, 88 
of circle, 90 

Vector product of two vectora, 
20 

non-commutative property, 
11 

distributive law, 24 
in rectangular co-ordinatca, 

24 

differentiation of, 31 
divergence of, 49 


Vector quantities, i 
invariance of, 10 
Vector triple product of vectorSj 
27 

Vectors, addition and subtrac¬ 
tion of, 6 

associative law for sum, 8 
components of, 8 
sum in rectangular co-ordiii' 
ates, 11 

products of t^vo, 14 
scalar, 14 
vector, 20 
triple products, 37 
differentiation of, 29 
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